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Preface 


The  idea  that  nitrogen  doping  strengthens  a  silicon 
lattice  is  intriguing  in  light  of  the  continuous  search  for 
better  semiconductor  materials.  A  review  of  the  literature 
showed  that  strengthening  effect  of  nitrogen  has  been 
studied,  but  the  electrical  parameters  had  not.  This  work 
is  an  effort  to  fill  that  gap. 

Dr.  W.  C.  Mitchel  of  the  Materials  Laboratory  ( AFWAL ) 
proposed  this  project.  He  supplied  the  material  and 
provided  the  Hall  effect  system  on  which  the  experiments 
were  performed.  I  am  grateful  for  his  suggestions  on  how 
to  proceed  with  the  experiments  and  his  advice  in  analyzing 
the  data.  I  would  like  to  thank  Lynnette  Brown  who  gave  me 
much  needed  help  with  the  day  to  day  operation  of  the 
equipment. 

I  would  like  to  thank  CAPT  D.  Kitchen,  my  advisor,  for 
his  guidance  and  encouragement  and  also  MAJ  W.  Sutton,  my 
reader,  for  his  suggestions  in  editing  this  work. 

Finally,  I  express  my  gratitude  to  my  husband,  Michael, 
for  his  patience  and  help  in  putting  together  the  final. 

Mary  A.  O'Leary 
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Abstract 


The  electrical  parameters  of  a  boule  of  nitrogen  doped 
float  zone  silicon  were  studied  by  Hall  effect  analysis. 
Samples  were  annealed  at  temperatures  ranging  from  800  to 
900  °C.  In  addition,  a  section  of  this  boule  was  neutron 
transmutated  (NTD)  to  increase  the  n-type  doping.  Samples 
from  the  NTD  section  were  annealed  at  800  ^C. 

Resistivity  and  mobility  varied  considerably  from 
sample  to  sample,  but  the  variation  is  not  a  function  of 
annealing  temperature.  The  annealed  Si:N  samples  were 
found  to  be  inhomogeneous;  however,  the  Si:N  NTD  samples 
were  homogeneous.  In  addition,  annealing  activated  deep 
energy  levels.  The  shallow  energy  levels  reported  in 
studies  on  ion  implanted  nitrogen  in  silicon  were  not 
found.  The  only  shallow  energy  level  found  was  phosphorus. 

The  conclusion  is  nitrogen  does  not  go  into 
substitutional  sites  in  silicon  to  any  great  extent. 
Although  what  happens  to  nitrogen  in  silicon  is  unknown  it 
is  suggested  that  nitrogen  may  form  silicon-nitride 
complexes  and  precipitates.  This  could  be  the  mechanism 
for  strengthening  the  lattice. 
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ELECTRICAL  PROPERTIES  OF  BULK  GROWN 
NITROGEN  DOPED  SILICON 


I.  Introduction 


Single  crystal  silicon  is  by  far  the  most  widely  used 

semiconductor  material  in  the  world  today.  The  dominance 

of  silicon  over  other  materials  is  mostly  due  to  its 

material  properties  but  probably  equally  important  is  the 

fact  that  high  quality  single  crystal  silicon  is  cheaper 

and  easier  to  grow  than  any  other  semiconductor  material. 

As  the  demand  for  higher  device  density  and  faster 

performance  increases,  manufacturers  are  pushed  to  produce 

silicon  with  fewer  lattice  defects.  The  control  of  defects 

in  a  semiconductor  substrate  is  important  since  defects 

directly  effect  device  characteristics  and  can  cause  device 

failure.  If  device  manufacturers  could  make  their  product 

as  small  as  they  liked  without  regard  to  the  limits  of 

lithography  and  other  processes: 

. .  localized  defect  clusters  in  an 

otherwise  well  controlled  silicon  wafer  may 
well  become  the  limiting  factor  for  improved 
device/circuit  performance  and  product 
yield.  This  may  especially  be  the  case  for 
VLSI  circuitry  as  the  probability  of  a 
defect  occurring  within  a  die  increases  as 
the  die  size  increases  (for  constant  defect 
density)  while  concurrently,  the  probability 
of  a  defect  resulting  in  the  occurrence  of  a 


circuit  malfunction  increases  as  the 
minimum  lithographic  dimension  decreases 
(16:90). 


Lithographic  techniques  are  improving,  making  smaller 
geometries  possible.  However,  a  clear,  sharp  image  relies 
on  a  flat  surface.  "Today's  technology  produces  a  flatness 
of  less  than  6  microns;  however,  flatness  as  low  as  2 
microns  or  less  will  be  required  by  1985,  and  less  than  1 
micron  by  1990"  (35:101).  Many  device  processing  steps 
require  high  temperatures  (up  to  1200  °C)  which  can 
result  in  wafer  warpage  due  to  residual  lattice  strain. 

Furthermore,  in  order  to  maximize  charge  carrier 
mobility,  which  is  the  parameter  that  ultimately  limits 
device  speed,  the  crystal  lattice  must  be  as  defect  free  as 
possible.  That  is,  the  crystal  needs  to  be  free  of 
unwanted  impurity  atoms  and  other  irregularities  which 
disrupt  the  periodicity  of  the  lattice,  such  as  dislocation 
defects.  The  dilemma  facing  very  large  scale  integrated 
circuit  (VLSI)  and  very  high  speed  integrated  circuit 
(VHSIC)  device  manufacturers  is  that  high  purity, 
dislocation  free  silicon  warps  badly  during  device 
processing  at  elevated  temperatures.  It  has  recently  been 
discovered  that  the  oxygen  that  contaminates  electronic 
grade  commercial  silicon  is  responsible  for  strengthening 
the  wafers  (30:L51).  Oxygen  contamination  is  a  consequence 
of  the  Czochralski  (CZ)  method  of  growing  a  single  silicon 
crystal.  Ninety  percent  of  the  silicon  substrates  in  use 
today  are  manufactured  by  the  CZ  method.  The  other  ten 
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percent  of  the  single  crystal  silicon  is  manufactured  by 
the  float  zone  (FZ)  method,  usually  for  applications 
requiring  high  purity  silicon.  Typically,  FZ  silicon  has 
1000  times  less  oxygen  than  CZ  silicon.  FZ  grown  silicon 
has  been  rejected  as  a  substrate  for  nearly  all  scales  of 
integration  because  of  its  poor  resistance  to  warping 
during  processing.  Although  the  oxygen  has  the  beneficial 
effect  of  strengthening  the  wafer  by  pinning  the 
dislocation  faults,  it  is  at  the  expense  of  creating 
stacking  faults  and  decreasing  carrier  mobilities. 

In  response  to  this,  new  ideas  are  being  explored  to 
either  control  or  exploit  crystal  defects.  The  intrinsic 
gettering  approach  exploits  the  oxygen  in  CZ  silicon  by 
precipitating  this  element  onto  defects  in  the  bulk  of  the 
wafer.  This  leaves  a  defect  free  zone  on  the  surface  for 
device  processing.  The  intrinsic  gettering  process  adds 
extra  processing  steps  to  device  fabrication.  It  is  a 
complex  process  requiring  several  high  temperature  cycles, 
and  results  in  a  wafer  with  altered  characteristics  (10). 

In  addition  to  this,  oxygen  is  difficult  to  control  during 
crystal  growth  (11). 

In  1 981  a  group  of  Japanese  researchers  doped  a  boule 
of  FZ  grown  silicon  with  nitrogen.  They  reported  that  this 
material  had  fewer  dislocations  defects  than  ordinary  FZ 
silicon  and  slightly  less  than  CZ  silicon  after  two  heat 
treatments  (1:63).  A  reproduction  of  the  X-ray  topographs 
of  the  wafers  is  shown  in  Figure  1. 
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Figure  1 .  Effect  of  Nitrogen  on  Dislocation 
Nucleation  and  Movement  in  Silicon 
(1:71 ) 


The  nitrogen  concentration  of  the  boule  was  1.5  x 


1015  cm"3.  This  is  1000  times  less  than  the  oxygen 
doped  CZ  wafer  at  1.7  x  101®  cm-3.  The  conversion 
efficiency  of  donors  by  nitrogen  is  2  to  5%,  thus  nitrogen 
is  essentially  electrically  inert.  This  suggests  that  the 
carrier  mobilities  could  be  higher  than  CZ  silicon  and  that 
nitrogen  doped  wafers  will  survive  thermal  processing. 
Nitrogen  doped  silicon  (Si:N)  would  then  be  suitable  for 
VLSI  and  VHSIC  applications. 

Approach 

While  the  mechanical  properties  of  FZ  Si:N  have  been 

reported  in  the  literature,  the  electrical  properties  have 

not.  Matsato  Imai  and  Koji  Sumino  reported  in  their  X-ray 

topographic  studies  of  dislocation  mobility  in  silicon  that 

"Nitrogen  atoms  are  concluded  to  have  a  strong 
immobilization  effect  for  dislocations  in  silicon 
crystal  since  the  cessation  of  motion  occurs  at  a 
concentration  as  low  as  5.4  x  10^5  cm"3  (13:619)." 

Photoluminescence  (PL)  studies  have  been  undertaken  by 

various  other  groups  (38,39,40).  These  studies  have 

produced  data  about  the  activation  energy  of  nitrogen  in 

silicon,  deep  level  trap  states,  and  the  conversion 

efficiency,  that  is,  the  ratio  of  electrically  active 

nitrogen  to  the  total  concentration  of  nitrogen  in  the 

material.  The  electrical  parameters  of  ion  implanted 

nitrogen  in  p-type  silicon  have  been  studied. 

Unfortunately,  there  are  disagreements  as  to  what  the 


activation  energy  is,  what  the  deep  level  activation  energy 
is,  and  whether  the  shallow  level  is  caused  by 
substitutional  nitrogen  or  silicon  nitride  complexes. 

In  this  study,  the  electrical  parameters  of  bulk  grown 
Si:N  were  studied  using  the  Hall  Effect  technique.  This 
technique  is  frequently  used  to  characterize  semiconductor 
materials.  The  Hall  Effect  technique  gives  a  wealth  of 
data  about  the  material  being  studied.  Resistivity, 
carrier  concentration,  and  carrier  mobility  are  determined 
from  room  temperature  data.  If  the  experiment  is  performed 
over  a  wide  range  of  temperatures,  such  as  liquid  helium 
temperature  (4  Kelvin)  to  room  temperature  (approximately 
300  Kelvin),  the  activation  energy  levels  and  doping 
concentration  can  be  determined. 

Experiments  Performed 

Samples  of  Si:N  were  annealed  at  various  temperatures 
from  800  to  900  °C  and  compared  an  unannealed  sample  and 
a  CZ  sample  annealed  at  900  °C.  In  addition,  a  two  inch 
chunk  of  Si:N  was  sent  to  the  University  of  Missouri 
Nuclear  Reactor  for  neutron  transmutation  doping  (NTD).This 
technique  increased  the  n-type  doping  of  the  Si:N  material 
by  the  nuclear  reaction  (22:1): 

30Si(n,V)31Si  - —  31P  +  0-  (1) 

where  n  »  neutron 
y  *  gamma 
0- =  beta 

"....and  the  fact  that  the  3^si  isotope,  with  an 
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abundance  of  3.1%,  is  uniformly  distributed  through¬ 
out  the  crystal  lattice  (17:139)". 

The  data  from  the  NTD  Si:N  was  compared  to  the  data  from  a 

previous  experiment  performed  by  W.C.  Mitchel  of  AFWAL/MLPO 

on  a  high  purity  FZ  sample  that  was  also  neutron  trans- 

mutated  (23). 

The  Hall  effect  experiments  were  performed  at  the 
Materials  Laboratory  (AFWAL)  on  the  Hall  effect  systems 
owned  by  the  Laser  and  Optical  Materials  Branch.  The 
resulting  data  was  analyzed  by  a  computer  program  written 
specifically  for  use  with  the  Hall  experiment.  The  program 
resides  on  the  Prime  computer  at  the  Materials  Laboratory. 

The  Hall  Effect  experiments  only  give  mobility  data  on 
the  majority  carrier.  Because  nitrogen  is  a  donor  in 
silicon,  the  mobilities  obtained  were  for  electrons  only. 

In  order  to  study  the  mobility  of  holes  the  material  would 
have  to  be  p-type.  It  is  possible  to  ion  implant  n-type 
substrate  with  a  p-type  ion.  By  reverse  biasing  the 
pn-j unction  formed  by  the  p-type  ion  implanted  layer  and 
the  n-type  substrate,  the  p-type  layer  is  electrically 
isolated  from  the  substrate  so  that  the  Hall  effect 
experiment  can  be  performed.  However,  in  order  to  study  an 
ion  implanted  layer,  the  system  has  to  be  modified  to 
include  and  extra  power  supply  and  ammeter.  A  detailed 
discussion  of  this  is  provided  in  Appendix  A. 


Results 


The  results  of  the  Hall  experiments  on  these  samples  of 
Si:N  showed  that  the  material  is  not  homogeneous.  The 
electrical  parameters  varied  widely  from  sample  to  sample. 
The  electrical  parameters  do  not  appear  to  change  as  a 
function  of  the  annealing  temperature.  Although,  annealing 
does  activate  deep  energy  levels  in  the  material.  In 
addition,  annealing  can  cause  a  reasonably  homogeneous 
sample  to  become  inhomogeneous.  The  NTD  Si:N  samples  were 
dominated  by  the  phosphorus  doping.  Annealing  did  not 
result  in  inhomogeneity.  The  energy  levels  reported  in  the 
previously  published  PL  studies  were  not  found  in  these 
experiments.  A  PL  study  of  samples  from  the  same  boule  of 
Si:N  used  in  the  Hall  experiments  did  confirm  the  presence 
of  nitrogen  (6). 

This  series  of  experiments  showed  that  the  electrical 
properties  of  Si:N  are  not  consistent  due  to  the 
inhomgeneous  nature  of  the  material.  Thus,  without  further 
study  on  how  to  overcome  this  problem  Si:N  cannot  be 
recommended  as  a  VLSI  or  VHSIC  substrate. 

Organization  of  Thesis  Chapters 

In  order  to  understand  how  these  conclusions  were 
reached  the  remaining  chapters  give  an  overview  of  silicon 
and  current  issues  in  manufacturing  an  explanation  of  the 
Hall  experiment,  and  a  detailed  discussion  of  the  results. 
Specifically  Chapter  II  has  a  discussion  on  the  silicon 
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crystal  lattice  and  defects ,  methods  of  manufacturing,  and 
previous  work  on  nitrogen  in  silicon.  Chapter  III  explains 
the  Hall  effect,  how  the  electrical  parameters  are  derived 
from  the  experimental  data,  and  an  explanation  of  the 
computer  analysis  of  the  data.  Chapter  IV  describes  the 
experimental  procedure  and  the  Hall  effect  system.  The 
results  of  the  experiments  are  given  in  Chapter  V,  while 
Chapter  VI  gives  an  analysis  of  the  results.  Chapter  VII, 
the  last  chapter,  gives  the  conclusions  and  recommend¬ 


ations 


II.  Background 


Since  crystal  defects  affect  device  perf  ounance ,  the 
various  types  of  defects  and  how  they  occur  need  to  be 
explained.  Certain  defects  occur  when  the  crystal  is 
grown,  others  occur  during  device  processing.  The  various 
types  of  defects  will  be  described  first,  then  the  methods 
of  manufacturing  single  crystal  silicon.  With  this 
background  the  implications  of  intrinsic  gettering  and 
nitrogen  doped  FZ  Si  can  be  understood.  From  this, 
previous  work  on  the  effect  of  nitrogen  in  silicon  is  more 
easily  explained. 

Lattice  Defects 

Solid  state  semiconductor  theory  is  based  on  taking 
advantage  of  defects  in  a  lattice.  Doping  atoms,  such  as 
phosphorus  or  boron,  add  extra  charge  carriers  to  the 
material.  Phosphorus  adds  extra  electrons  (negative  charge 
carriers),  while  boron  atoms  create  a  deficiency  of 
electrons,  called  holes  (positively  charge  carriers).  When 
a  material  is  doped  to  have  excess  negative  carriers 
(n-doped)  the  electrons  are  referred  to  as  majority 
carriers  while  the  holes  are  minority  carriers.  For 
p-doped  material  (excess  holes)  the  holes  are  the  majority 
carriers  while  the  electrons  are  the  minority  carriers. 

The  electrical  characteristics  of  devices  are  described  by 


the  behavior  of  the  carriers.  The  flow  of  the  charged 
carriers  through  the  material  is  an  electrical  current. 
When  an  electron  recombines  with  a  hole,  the  current 


decreases.  Unfortunately,  unintentional  impurities,  such 
as  oxygen,  and  lattice  deformations  provide  extra 
generation/recombination  sites.  Since  they  are  unwanted 
they  are  called  a  traps.  Traps  reduce  minority  carrier 
lifetime  and  diffusion  length.  The  equations  governing 
minority  carrier  lifetime  and  diffusion  length  are 
(8:124,1 34) : 


where 


tm  =  1/(s  x  vth  x  Nt)  (seconds)  (2) 

Lm  =  (D  x  tm) _1 /2  (cm)  ( 3 ) 

tm  =  minority  carrier  lifetime  (sec) 
vth  =  thermal  velocity  =  (3KT/M)1/2  (cm/sec) 
s  =  effective  trap  cross  section  (cm2) 

Nt  =  concentration  of  traps  (cm~3) 

Lm  =  diffusion  length  of  minority  carriers  (cm) 

D  =  diffusion  coefficient  of  minority  carriers 
(cm2/ sec) 

Mobility  is  mainly  dependent  on  charged  center  scattering. 
Traps  are  charged,  so  they  scatter  electrons  which  reduces 
mobility  (29:97-102,253-263).  Defects  are  not  uniformly 
distributed;  therefore,  minority  carrier  lifetime  is  not 
constant  within  the  wafer  (Figure  2).  This  lack  of 
homogeneity  can  result  in  differences  in  the  the  electrical 
parameters  of  the  individual  devices  in  an  integrated 
circuit  chip. 

Defects  occur  as  the  crystal  is  grown  or  during  device 
processing.  Grown-in,  or  intrinsic,  defects  include  the 
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various  types  of  point  defects,  dislocations,  and  stacking 


faults.  Process  induced,  or  extrinsic,  defects  include 
surface  impurities,  dielectric  definition,  design  related 
faults,  and  process  generated  dislocation  and  stacking 
faults.  Extrinsic  defects  are  controlled  by  following 
procedures  for  cleanliness  and  using  proper  design  rules. 
Issues  in  dielectric  definition  are  beyond  the  scope  of 
this  paper.  Process  induced  dislocation  and  stacking 
faults  are  directly  related  to  the  quality  of  the  starting 
material . 

Defects  disrupt  the  periodicity  of  the  lattice  and 
produce  stresses  in  the  immediate  vicinity,  producing 
traps  and  scattering  centers.  Point  defects  disrupt  one 
lattice  site,  but  dislocations  and  stacking  faults  disrupt 
whole  sections. 


Figure  2.  Minority  Carrier  Lifetime  Distribution  (28:331) 
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The  diamond  structure  of  crystalline  silicon  is  rather 
open,  so  impurity  atoms  can  occur  interstitial ly  as  well  as 
subst i tut iona 1 ly .  Figure  3  gives  examples  of  each  of  the 
point  defects.  Note  that  silicon  atoms  can  also  occur 
between  lattice  sites,  which  is  referred  to  as  a 
self -interstitial .  The  oxygen  in  silicon  is  primarily 
interstitial.  When  an  atom  of  silicon  is  missing  from  the 
lattice  its  called  a  Schottky  defect  or  a  vacancy,  but  if 
an  atom  leaves  its  lattice  site  and  becomes  interstitial 
the  combination  of  a  vacancy  and  an  interstitial  is  called 
a  Frenkel  defect. 

Dislocation  and  stacking  faults  are  extra  planes  of 
atoms  in  the  lattice.  Figure  4  illustrates  an  edge 
dislocation.  The  lattice  is  in  compression  above  the 
the  edge  and  in  tension  below  it.  A  screw  dislocation 
occurs  when  an  entire  block  of  atom  slips  out  of  place  as 
shown  in  Figure  5.  This  produces  a  twisting  stress  on  the 
lattice.  A  stacking  fault.  Figure  6,  can  be  any 
combination  of  edge  and  screw  faults. 

Point  defects  and  dislocations  can  slip  and  -ove  within 
the  lattice,  especially  when  heated.  In  fact,  Frenkel  and 
Schottky  defects  have  an  activation  energy.  Dislocation 
faults  are  particularly  troublesome.  If  an  epitaxial  layer 
is  grown  over  a  dislocation  fault,  the  dislocation  wiJl 
grow  into  the  new  layer.  Dislocations  can  be  generated 
during  high  temperature  processing  and  will  move  and  grow 
when  stressed. 
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Figure  3.  Schematic  Illustration  of  Point  Defects: 

Substitutional  Boron,  Substitutional  Carbon 
Interstitial  Oxygen,  Silicon  Self-Inter¬ 
stitial,  Vacancy  or  Schottky,  Substitu¬ 
tional  and  Interstitial  Copper,  and  Frenk¬ 
el.  Note  the  surface  steps  of  silicon 
atoms .  (11:212) 


Figure  4.  Edge  Dislocation  (23:212) 


Figure  5.  Screw  Dislocation  (16:154) 


Figure  6.  Stacking  Fault 


The  dislocation  density  on  the  surface  of  a  wafer  is 
revealed  by  X-ray  topography.  Dislocation  areas  are  more 
susceptible  to  etching  because  the  region  is  stressed. 
Illuminating  the  surface  with  X-rays  shows  the  dislocations 
as  spots.  Typically,  electronic  grade  silicon  has  less 
than  100  dislocations/cm^ ,  although  less  than  1 O/cm^  is 
possible  (23).  Thus,  most  single  crystal  silicon  is  called 
dislocation  free.  Even  so,  any  dislocations  are  becoming 
worrisome  to  device  manufacturers,  because  of  ever 
decreasing  component  size  and  the  simultaneous  increase  in 
chip  size. 

Single  Crystal  Manufacturing 

Commercially  available  single  crystal  silicon  is  made 
by  the  Czochralski  Puller  method  or  the  float  zone  method. 
Although  FZ  crystals  are  purer  than  CZ  crystals,  CZ 
crystals  account  for  90%  of  the  world's  production,  while 
FZ  is  10%. 

In  the  CZ  method  polycrystalline  silicon  is  melted  in  a 
quartz  crucible.  The  crucible  is  placed  in  a  graphite 
susceptor.  A  small  seed  crystal  of  the  proper  crystal 
orientation  is  introduced  to  the  surface  of  the  melt,  then 
is  slowly  pulled  out.  As  the  seed  crystal  is  pulled  out  of 
the  melt,  new  material  is  deposited  which  continues  the 
same  lattice  configuration.  Figure  7  shows  this 
arrangement  schematically.  The  pull  rate  is  in  the  range 
of  2  to  5  x  1 0“^3  cm/sec.  The  diameter  of  the  boule  is 
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determined  by  the  temperature  and  the  pull  rate.  The 
largest  boule  reported  is  12  inches  in  diameter  (23). 

Besides  the  doping  atoms  which  are  intentionally  added, 
the  melt  is  contaminated  by  oxygen  and  carbon  atoms.  The 
hot  melt  slowly  dissolves  the  crucible,  then  breaks  the 
silicon  dioxide  down  into  silicon  and  oxygen.  Carbon  atoms 
from  the  susceptor  contaminate  the  atmosphere  surrounding 
the  melt.  A  few  of  these  atoms  dissolve  into  the  melt. 
Carbon  contamination  varies  between  0.1  to  3  parts  per 
million  atoms  (ppma)  or  equivalently  0.05  to  1 . 5  x  10^ 


Figure  7.  Schematic  of  Czochralski  Puller  Method  (11:92) 
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atoms/cm  (34:113).  The  role  of  carbon  Is  uncertain,  but 
seems  to  aid  in  the  precipitation  of  oxygen. 

Typically  oxygen  concentration  in  Cz  crystal  varies 
from  18  to  40  ppma  (0.9  to  2  x  10^  atoms/cm3)  along 
the  length  of  the  boule  (34:112-113).  Efforts  to  control 
the  oxygen  concentration  can  reduce  axial  variation  to 
between  24  to  38  ppma  (1.2  to  2  x  1  01  8  atoms/cm3). 

Applying  a  magnetic  field  across  the  melt  to  dampen  out 
convection  currents  reduces  the  oxygen  variation  and 
concentration  to  10  to  14  ppma  (5  to  7  x  1017 
atoms/cm3).  Figure  8  shows  the  oxygen  distribution  from 
top  to  bottom  of  the  boule  for  three  different  types  of  CZ 
silicon.  The  distribution  of  oxygen  across  a  wafer  is  also 


Figure  8.  Axial  Distribution  Figure  9.  Resistivity 

of  Oxygen  in  Three  CZ  Process-  (Dotted)  and  Oxygen 
es  (34:112)  Concentration  Distribu¬ 

tion  (Solid)  (1:67) 
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not  constant.  Figure  9  shows  the  oxygen  distribution  and 
resistivity  across  the  surface  of  three  different  wafers. 
Single  crystal  silicon  of  very  high  purity  is  made  by 


the  Float  Zone  method.  Most  FZ  silicon  is  made  for 
applications  where  high  resistivity,  9  to  30  K  ohms-cm,  is 
necessary,  such  as  infrared  detectors.  Examples  where  very 
pure  crystals,  but  not  necessarily  high  resistivity,  are 
needed  are  power  devices  and  devices  such  as  discrete 
transistors  and  diodes.  Most  FZ  silicon  is  grown  for  these 
applications . 

The  polycrystalline  boule  is  held  vertically  as 
schematically  illustrated  in  Figure  10.  An  RF  heating 


source  heats  the  bottom  while  a  seed  crystal  is  brought  up 
from  underneath.  The  heat  source  is  slowly  moved  upward 


Figure  10.  Schematic  of  Figure  11.  Float  Zone  with 
Float  Process  (21:376)  Needle-Eye  Coil  (17:138) 


and,  as  it  goes,  the  polycrystalline  material  takes  on  the 
lattice  structure  of  the  seed  crystal.  The  surface  tension 
of  the  liquid  keeps  the  melted  silicon  from  flowing  away. 
Figure  11  is  a  schematic  the  current  FZ  process  with  a 
small  "needle-eye"  RF  induction  coil. 

Since  there  is  no  crucible  to  contaminate  the  melt,  the 
oxygen  content  of  FZ  silicon  can  be  as  low  as  0.01  to  0.06 
ppma  (0.5  to  3  x  10^  atoms/cm^).  Because  there  is  no 
graphite  susceptor,  the  carbon  concentration  in  FZ  silicon 
is  less  than  0.01  ppma  (5  x  10^  atoms/cm^). 

Controlling  Defects 

Controlling  defects  is  of  great  concern  to  the 
semiconductor  industry.  The  role  of  oxygen  has  been 
studied  extensively  as  it  is  a  major  contaminant.  Until 
recently,  oxygen  was  thought  of  as  an  unfortunate 
consequence  of  the  CZ  process.  Efforts  have  been  made  to 
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reduce  the  oxygen  concentration  in  the  CZ  silicon;  however, 
it  is  now  realized  that  oxygen  hardens  the  wafer  by  pinning 
dislocation  faults  and  keeping  them  from  multiplying. 

The  silicon  industry  in  now  making  efforts  to  exploit 
the  benefits  of  oxygen.  Although  oxygen  pins  dislocation 
faults,  it  does  so  by  forming  stacking  faults.  Thus, 
faults  are  not  eliminated  but  controlled.  If  the 
concentration  of  oxygen  is  greater  than  30  ppma  (1.9  x 
1018  atoms/cm8),  SiyOx  precipitates  form,  leading 
to  wafer  warpage  during  thermal  processing. 

If  a  defect  free  zone  is  produced  on  the  surface  of  a 
wafer  and  active  electronic  devices  are  confined  to  this 
region,  then  it  would  not  matter  if  the  bulk  has  defects. 
This  new  idea  is  called  intrinsic  gettering.  Basically  it 
consists  of  three  temperature  cycles  ranging  from  650  to 
1200  °C  for  6  to  36  hours.  This  produces  a  defect  free 
or  denuded  layer  10  to  40  microns  deep  (Figure  12). 

However,  care  must  be  taken  when  planning  the  thermal 
cycles  for  intrinsic  gettering.  D.  Huber  and  J.  Reffle 
caution  that  instrinsic  gettering  introduces  a  new 
parameter  to  device  processing  (10:139).  They  carried  out 
experiments  which  show  that  if  the  oxygen  has  precipitated, 
the  wafer  is  subject  to  warping  during  subsequent 
processing  (Figure  13).  As  an  aid  to  designing  the 
appropriate  temperature  cycles  for  the  intrinsic  gettering 
process,  Huber  and  Reffle  processed  a  group  of  wafers  at 
1100  °C  for  10  hours  then,  exposed  to  different 
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Figure  13.  Warpage  Due  to  Stress  Cycles  at  1050  °C 
(10:140) 

nucleation  (pre-anneal ) and  precipitation/gettering 
(out-diffusion)  cycles.  Subsequently,  bipolar  transistors 
were  made.  "After  initial  oxidation  no  defects  are 
delineated;  however,  minority  carrier  diffusion  length 
measurements  show  a  drop  from  the  initial  values  of  250  to 
500  microns  down  to  values  between  50  and  100  microns" 
(10:141).  The  beneficial  effects  of  intrinsic  gettering 
may  be  negated  if  wafers  are  not  properly  processed. 


The  group  of  Japanese  researchers  (T.  Abe  et  al.)  that 
reported  that  nitrogen  also  pins  dislocation  faults  are 
continuing  their  investigations  (2,13,31,32,33).  When  it 
was  discovered  that  oxygen,  a  contaminant,  strengthened  C7, 
silicon  wafers,  the  Japanese  investigated  other  dopants 
that  have  the  same  effect.  Matsato  Imai  and  Ko j i  Sumino 
reported  that  oxygen,  nitrogen,  and  phosphorus  have  a 
locking  effect,  but  carbon  does  not  (13:618-619).  Imai  and 
Sumino  used  in-situ  x-ray  topography  to  study  the 
propagation  and  movement  of  the  dislocations  on  the  surface 
of  silicon. 

Apparently  this  group  of  researchers  are  looking  for  a 
material  that  has  a  strong  locking  effect  on  silicon  that 
is  either  electrically  inert  or  doped  in  such  low 
concentration  that  it  does  not  affect  the  silicon.  This  is 
a  new  area  of  investigation.  Literature  on  the 
strengthening  effect  of  nitrogen  starts  around  1980. 
Although  nitrogen  does  seem  to  have  good  properties  in 
silicon,  it  is  primarily  used  for  discrete  devices  and 
power  devices  (26:66). 

Previous  Nitrogen  Studies 

All  group  V  elements  in  the  periodic  table  are  donors 
in  silicon  (group  IV);  however,  the  behavior  of  nitrogen  is 
the  least  understood.  Early  attempts  to  dope  silicon  with 
nitrogen  using  conventional  doping  techniques  were  un¬ 
successful.  Nitrogen  reacts  with  molten  silicon  to  give 


silicon  nitride  (S13N4)  which  readily  dissolves  into 

the  melt  (  1  5 :  4  3  1  )  .  In  1973,  Yutsurgui  detei  mined  the  phase 

diagram  of  N-Si  in  extremely  low  nitrogen  concentrations. 

!!'•  discovered  tha*  a  small  amount  of  nitrogen  can  ho 
present  in  silicon  and  reported  that  the  solubility  limit, 
of  nitrogen  in  silicon  is  4.5  +_  1.0  x  10^  cm~3  in 
solid  silicon  and  about  6  x  10^®  cni-^  in  liquid  silicon 
(43:978),  making  nitrogen  the  least  soluble  group  V  element 
in  silicon.  Further,  from  his  infrared  absorption  studies 
he  reported  that  less  than  1  %  of  the  nitrogen  goes  into  the 
lattice  substitutionally ,  and  concluded  that  most  of  the 
nitrogen  remains  un-ionized  (43:978). 

The  difficulties  of  introducing  nitrogen  into  silicon 
was  well  known  before  Yatsurgai's  results  were  published. 
Other  groups  got  around  this  difficulty  by  ion  implanting 
nitrogen  into  p-type  silicon  substrate.  One  of  the 
earliest  studies  of  ion  implanted  nitrogen  was  performed  in 
1968  by  Zorin  et  al.  of  the  USSR  (44).  They  were  able  to 
isolate  the  n-type  nitrogen  layer  by  reverse  biasing  the 
pn-junction  formed  by  the  n-type  ion  implanted  layer  and 
the  p-type  substrate,  and  performed  the  Hall  experiment  on 
it.  They  reported  the  activation  energy  of  nitrogen  to  be 
0.045  eV  below  the  conduction  band,  that  ion  implanted 
nitrogen  occupies  substitutional  sites,  and  the  nitrogen 
behaves  like  a  normal  donor  (44:111,112).  This  result  is 
remarkable  in  that  the  activation  energy  of  phosphorus  is 
also  0.045  eV  below  the  conduction  band,  and  the  donor 
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concentration  measured  was  less  than  1%  of  the  implanted 
nitrogen  ions. 

A  similar  Hall  experiment  was  repeated  on  ion  implanted 

nitrogen  in  1975  by  J.  B.  Mitchell  et  al.  (25).  In  this 

experiment  FZ  and  CZ  samples  were  ion  implanted.  Mitchell 

reported  an  activation  energy  of  0.017  +.  0.002  eV.  In  this 

study  the  temperature  range  of  the  Hall  experiment  was  much 

greater,  from  20  K  to  300  K,  while  Zorin's  data  was  taken 

from  125  K  to  250  K.  From  a  previous  paper  of  atom-site 

location  measurements  (24),  Mitchell  found  that  greater 

than  95%  of  the  implanted  nitrogen  was  located  on 

interstitial  sites,  leaving  an  undetectable  amount  on 

possible  substitutional  sites. 

The  existence  of  a  unique  ionization  energy  which 
is  independent  of  anneal  time  and  implantation  temp¬ 
erature  (and  consequently,  the  initial  amount  of 
damage)  suggests  a  substitutional  position  for  the 
donor-active  nitrogen.  The  shallowness  of  the  level 
(0.17  +_  0.002  eV)  below  the  conduction  band  also 
suggests  a  substitutional  site,  as  all  known  defect 
centers  in  silicon  give  rise  to  deeper  states  in  the 
band  gap  (25:341). 

However,  he  also  suggests  that  there  may  be  a  defect 
conduction  mechanism  since  the  data  does  not  fit  the  usual 
curve  for  a  single  donor  level.  This  indicates  that  there 
may  be  some  form  of  defect  centers. 

In  a  later  study  (1980),  Keith  Brower  ion  implanted 
nitrogen  into  intrinsic  FZ  silicon  (4).  The  implanted 
layer  was  laser  annealed  to  restore  the  crystal  lattice. 
Brower  used  electron  paramagnetic  resonance  (EPR)  to  study 
the  structure  of  nitrogen  in  silicon.  He  concluded  that 


substitutional  nitrogen  and  neighboring  silicon  atoms  are 
" Jnhn  T<  •  1  1  < * r  d  i  st  or  I  <  •<]  off  I  In  •  i  i  I  <  •(  rahedi  . i  1  suhst  i  I  u I  i  t «i i . i  1 
sites  along  directions  which  preserve  axial  symmetry  about 
the  '  1  1  1  >  axis  (4:1628)."  Associated  with  the  distorted 
substitutional  nitrogen  is  an  activation  energy  oi  0.88  cV 
below  the  conduction  band.  Thermal  annealing  above  430 
°C  appears  to  modify  the  structure  of  the  distortion,  and 
"might  be  the  beginning  of  the  formation  of  silicon  nitride 
complexes  with  shallow  donor  levels  which  have  been 
observed  in  electrical  measu t ernents  by  others  (4:162CM." 

Since  the  discovery  of  the  strengthen ing  effect  of 
nitrogen  in  silicon,  other  studies  have  been  undertaken  to 
characterize  nitrogen  in  bulk  grown  silicon  rather  than  ion 
implanted  layers.  The  nitrogen  is  introduced  into  the 
boule  during  the  growth  phase  either  by  adding  nitrogen  gas 
to  the  argon  ambient  or  adding  silicon  nitride  to  the 
melt.  A  photoluminescence  (PL)  study  by  Tajima  et  al.  in 
1981,  revealed  the  presence  of  a  PL  line  at  1.1223  +  0.0001 
eV  which  increases  in  amplitude  with  nitrogen  doping  (38: 
L423).  This  same  PL  line  was  also  found  in  ion  implanted 
nitrogen,  which  is  further  evidence  that  this  defect  is 
associated  with  nitrogen  (3:426).  PL  can  be  used  to 
determine  the  impurity  concentration  of  shallow  impurities 
in  bulk  material.  The  nitrogen  concentration  in  this  study 
was  determined  by  activation  analysis.  Tajima  considers 
the  1.1223  eV  line  to  be  due  to  the  radiative  recombination 
of  an  exciton  bound  to  a  deep  nitrogen  center  (38:L425). 
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crystal,  but  not  necessarily  substitutional  nitrogen.  The 
discovery  of  the  1.1223  eV  line  defines  a  non-destructive 
test  tor  the  presence  of  nitrogen  in  a  sample. 

Another  experiment,  performed  a  year  later  by  Tokumaru 
et  al. ,  studied  the  material  using  deep  level  transient 
spectroscopy  ( DLTS )  in  an  eflort  to  find  the  0.58  eV  deep 
level  reported  by  Brower  (40).  Instead  of  the  0.58  eV 
level,  two  other  deep  Levels  were  found:  0.19  eV  at  0.1 
of  the  nitrogen  doping  concentration,  and  0.28  eV  at  0.01% 
of  the  nitrogen  concentration.  However,  the  material 
Brower  used  was  ion  implanted  and  annealed  at  about  480 
°C.  The  nitrogen  doped  F 7,  crystals  used  in  this  study 
were  grown  at  a  much  higher  temperature,  thus  isolated 
substitutional  nitrogen  atoms  may  not  exist  in  this 
material  and  that  the  energy  levels  may  be  due  to  different 
types  of  nitrogen  related  complexes  (40:L444).  The  shallow 
donor  of  0.017  eV  may  exist  in  the  crystal,  but  is  too 
shallow  to  be  measured  by  DLTS. 

In  summary,  nitrogen  doped  silicon  has  been  studied  a 
number  of  different  ways,  each  reporting  a  different 
activation  energy  and  theory  for  nitrogen's  effect  on 
silicon.  Although,  there  is  evidence  for  substitutional 
nitrogen  in  ion  implanted  silicon,  there  is  no  clear  cut 
evidence  for  it  in  bulk  grown  silicon  or  that  the 
activation  energies  are  the  same.  Hall  effect  analysis  of 
this  material  has  not  been  published,  which  is  the  most 


sensitive  method  of  determining  electrical  parameters  in  a 
semiconductor.  The  experiments  undertaken  in  this  thesis 
used  the  previous  work  as  a  starting  point  to  determine  the 
activation  energies.  The  results  of  previous  studies  are 
summarized  in  Table  I. 


Table  I 

Comparison  of  Activation  Energies  and  Techniques  Used 


Researcher 

Year 

Material 

Energy  Level 

Method/Limitat ion 

Zorin 
et  al. 

1  968 

Ion 

Implanted 

0.045  eV 

Hall:  temperature 

range  100  -250  K 
too  limited  to  find 
much  shallower  or 
deeper  levels 

Mitchell 
et  al. 

1  974 

Ion 

Implanted 

0.017  eV 

Hall:  temperature 
range  20-250  K, 
not  high  enough  to 
find  deep  levels 

Brower 

1980 

Ion 

Implanted 

0.58  eV 

EPM:  cannot  detect 
shallow  levels 

Tokumaru 
et  al. 

1982 

Bulk 

0.19  eV 

0.28  eV 

DLTS :  cannot  detect 
shallow  levels 

III.  Tho  ila  1 1  Effect 


The  Hall  effect  experiments  were  performed  at  the  Laser 
and  Optical  Materials  Branch,  Materials  Laboratory  on 
equipment  developed  by  P.  Hemenger  specifically  for  testing 
semiconductor  material  (9).  However,  before  describing  the 
equipment  in  detail,  the  Hall  effect  will  be  explained. 

This  will  give  the  reader  an  understanding  of  the 
principles  involved  in  the  experiment  and  the  wealth  of 
information  that  results.  After  the  principles  are 
explained  the  computer  program  used  to  resolve  the  raw  dat> 
will  be  discussed. 

The  Hall  Effect 

When  a  magnetic  field  is  applied  to  moving  charges,  a 
force  is  induced  on  the  charges  that  is  mutually 
perpendicular  to  the  magnetic  field  and  the  direction  of 
charge  travel.  Thus,  charge  is  deflected  toward  the  edges 


Figure  14.  The  Hall  Effect 


of  the  material  with  opposite  polarities  building  up  on 
opposite  sides  (Fig.  14).  The  charge  builds  up  until  an 
electric  field  is  generated  that  balances  the  deflecting 
force  (27:2).  (This  phenomenon  was  discovered  in  18/9  by 
F.  H.  Hall).  The  voltage  induced  between  the  sides  of  the 
material  is  called  the  Hall  Voltage  ( ) . 

For  a  constant  magnetic  field  and  constant  current,  the 
Hall  Coefficient  is  defined  by: 


VH  =  Rj[  1 18/ 1  (volts)  |  4  ) 

where  I  =  current  (amps) 

B  =  magnetic  field  (gauss) 
t  =  thickness  of  sample  (cm) 

Rg  =  Hall  coefficient  (cm^/coul)  (7:12) 

The  Hall  coefficient,  RH,  is  a  material  parameter  that  is 

related  to  the  carrier  concentrations  and  the  mobilities. 

It  is  computed  directly  from  Equation  4: 

rH  =  Vjjt  /  ( IB )  (cm3/coul)  (5) 

when  conduction  is  dominated  by  one  type  of  charge  carrier, 

Rh  is  given  by: 

rH  =  -uHn/(enun)  for  electron  conduction  (6) 

or 

rH  *  UHp/(epup)  for  hole  conduction  (7) 

where 

uHn  =  Hall  mobility  of  electrons  (cm2/V-sec) 
uHp  =  Hall  mobility  of  holes  (cm2/v-sec) 
un  =  conductivity  mobility  of  electrons  (cm2/V-sec) 

Up  =  conductivity  mobility  of  holes  (cm2/V-sec) 
n  =  concentration  of  electrons  (cm“^) 
p  =  concentration  of  holes  (cm-*) 
e  =  charge  of  an  elect  on  (coul)  (19:2) 

The  ratio  of  the  Hall  mobility  to  the  conductivity 

mobility,  unn/un,  is  called  the  Hall  scattering  factor, 
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re  or  r-factor,  for  electrons.  (A  similar  scattering 
factor  is  defined  for  holes,  rp. )  The  scattering  factor 
is  a  function  of  the  band  structure  of  the  material, 
temperature,  magnetic  field  strength,  and  the  scattering 
mechanisms  (7:12-14).  All  material  related  effects  are 
contained  in  re.  Usually,  the  r-factor  can  only  be 
measured  over  a  limited  range.  It  can  be  calculated  or 
estimated,  but  in  most  cases  it  is  close  to  unity  and  this 
is  often  the  value  used  in  analysis.  In  this  work  the 
r-factor  is  assumed  to  be  unity,  thus,  Rjj  for  n-type 
material  is  reduced  to: 

Rh  =  r/en  =  1 /en  (cm3/coul)  (8 

By  re-arranging  Equation  8,  the  carrier  concentration  is 
simply: 

n  =  r/eRH  =  1 /eRH  (cm-3)  (9 

The  resistivity  is  determined  directly  from  measured 
voltage  current  characteristics  of  the  material.  Again, 
when  one  type  of  conduction  dominates,  resistivity  is 
related  to  the  carrier  concentration  and  conductivity 
mobility  by: 

pr  =  1 /enun  (ohm-cm)  (10 

where  pr  is  defined  as  resistivity  in  this  work.  Only 
%  and  pr  are  directly  measured.  The  Hall  electron 
mobility  is  defined  by: 

uHn  =  RH/Pr  (cm2/V-sec)  (11 

If  the  scattering  factor,  that  is,  the  r-factor  can  be 


determined  or  assumed  to  be  unity,  as  it  is  in  this  work, 
then  the  conductivity  mobility  is  given  by: 

un  =  uHn/r  (cm2/V-sec)  (12) 

or 

un  =  uHn  (cm2/V-sec)  (13) 

If  the  Hall  coefficient,  RH,  and  resistivity,  pr  of 
a  material  is  measured  over  a  wide  temperature  range,  then 
the  carrier  concentration  data  derived  from  these 
measurements  can  be  used  to  determine  the  electrically 
active  impurities  and  their  concentrations.  At  very  low 
temperatures,  for  example  near  20  K,  there  is  not  enough 
thermal  energy  to  ionize  the  impurities.  Thus,  carrier 
concentration  is  negligible  and  resistivity  is  high.  As 
the  temperature  is  increased  impurities  start  to  ionize, 
giving  up  carriers  to  the  conduction  or  valence  band. 
Finally  when  all  of  the  impurities  of  an  atomic  species 
become  ionized,  an  increase  in  temperature  does  not 
increase  the  carrier  concentration  and  the  sample  is  said 
to  be  saturated.  In  silicon  most  shallow  dopants  are 
completely  ionized  at  room  temperature  (300  K).  The  rate 
at  which  a  species  ionizes  with  temperature  is  related  to 
the  activation  energy.  The  activation  energy  is  determined 
from  the  slope  of  the  logarithm  of  the  carrier 
concentration  versus  inverse  temperature.  The  equation  for 
this  curve  is  derived  from  the  charge  balance  equation  and 
Fermi-Dirac  statistics  (19:12-17).  This  equation  contains 
the  species  concentration  and  the  activation  energy.  A 
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computer  program  can  use  this  complex  expression  to 
generate  a  least  squares  fit  curve  and  determine  the 
concentrations  and  activation  energies  for  a  number  of 
levels.  Thus,  Hall  effect  experiments  can  be  used  to 
determine  the  impurities  and  the  concentrations  in  a 
semiconductor . 


The  Curve  Fitting  Function  to  Determine  Carrier 
Concentration  and  Activation  Energy 


The  fitting  function  used  to  generate  a  least  squares 
fit  curve  to  the  carrier  concentration  points  versus 
inverse  temperature  is  derived  from  the  charge  balance 
equation : 

n  +]CNAk  =  P  +/LNDi  O' 

k  i 

where 

n  =  number  of  electrons  in  the  conduction  band 
p  =  number  of  holes  in  the  valence  band 
NDi  =  number  of  ionized  donors  (positively  charge.!) 

=  number  of  ionized  acceptors  (negatively 
charged)  (19:12) 

Sums  are  over  all  donor  and  acceptor  species.  If  the 

number  of  donors  greatly  exceeds  acceptors,  which  is  the 

case  for  n-type  silicon,  then  the  number  of  holes  in  the 

valence  band  is  negligible  and  equation  14  reduces  to: 

n  +LNAk  =5ZnD1  (IE 

k  i 

The  sum  over  the  acceptors  can  be  replaced  by  the  total 
number  of  acceptors,  N^,  because  the  number  of  holes  is 
negligible  at  room  temperature  and  below  and  all  of  the 
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acceptors  end  up  negatively  charged  (19:12).  If 
Fermi-Dirac  statistics  are  applied  to  the  donors  and 
Boltzmann  statistics  are  applied  to  the  electrons  in  the 
conduction  band  as  discussed  by  Putley  (27:122-137)  and 
Lang  and  Hemenger  (19:13-18,50-54),  the  charge  balance 
equation  becomes: 


where 

Nr  =  density  of  states  in  the  conduction  band 
g^j  =  degeneracy  of  the  excited  energy  states 
^Aij"  activation  or  ionization  energy  (eV) 
k  =  Boltzmann's  constant  (eV/Kelvin) 

T  =  temperature  (Kelvin)  ( 1 9 : 1 3 ; 27 : 1 32 ) 

The  density  of  states,  Nc,  is  roughly  dependent  on 
temperature  raised  to  3/2  power  ( 1 9: 1 7; 20: 4646 ) .  Unless 
the  sample  is  very  heavily  doped  or  the  activation  energy 
is  very  small,  the  excited  states  contribute  very  little  to 
the  sum  in  Equation  16,  so  that  becomes: 

Qi  =  [Nc  exp(-EAi/kT) ]/g  (18) 

where 

g  =  degeneracy  of  the  ground  state  (19:14) 

For  donors  in  silicon,  g  is  taken  to  be  two  (27:123;18). 

By  fitting  Equation  1 6  to  the  experimentally  measured 
temperature  dependent  carrier  concentration,  n(T),  the 
total  concentration  of  all  the  acceptors,  concentration  of 
each  donor,  and  the  activation  energy  of  each  donor  are 
determined  variables  of  the  fit.  As  an  example,  a  typical 
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Hall  experiment  performed  on  a  phosphorus  doped  silicon  is 
shown  in  Figure  15.  At  high  temperatures  all  of  the  donors 
are  ionized  thus,  n  =  Nq  -  NA.  In  the  high  temperature 
"exhaustion"  (or  saturation)  region  the  carrier 
concentration  is  independent  of  temperature  as  there  are  no 
more  electrons  to  be  given  up  by  the  donors.  This  can  be 
seen  at  the  left  side  of  Figure  15.  At  the  lowest 
temperatures  none  of  the  donors  are  ionized.  This  is  seen 
at  the  right  side  of  Figure  15.  The  electron  concentration 
is  given  by: 

n  =  (ND  -  NA)Q/NA  ( cm-3 )  (19) 

but  because  Q  is  exponentially  dependent  on  temperature,  n 
is  proportional  to  the  exponent  of  -EA/kT  (19:16).  By 
taking  the  logarithm  of  equation  19a  linear  equation 
results.  Plotting  the  logarithm  of  the  carrier  concen¬ 
tration  versus  the  inverse  of  the  temperature  (multiplied 
by  1000  as  shown  in  Fig.  15)  as  the  dependent  variable 
gives  the  activation  energy  is  the  slope  of  the  graph.  In 
the  example,  the  donor  concentration  is  3.19  x  lO1^ 
cm~3#  the  activation  energy  is  0.04445  eV,  and  the  number 
of  acceptors  is  1.33  x  10^2  cm--*.  The  activation 
energy  of  phosphorus  given  by  Sze  (37:21)  is  0.045  eV  which 
is  close  enough  to  identify  phosphorus  as  the  donor 
species. 

If  there  is  more  than  one  donor  species,  the  shape  of 
the  curve  is  more  complex  and  is  more  difficult  to 
resolve.  The  plot  of  carrier  concentration  versus  inverse 
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Figure  15.  Phosphorus  Doped  CZ  Silicon 
Annealed  at  900°C 


temperature  wiLl  still  liave  an  exhaustion  region  on  t!ie 
left  if  the  temperature  ranqe  of  the  experiment  was  high 
enough  to  ionize  all  donor;:,  but  the  analysis  of  the  lower 
temperature  range  will  depend  on  the  number  of  acceptors. 
The  shallowest  donors  will  not  be  seen  in  the  data  if  they 
are  overcompensated  by  the  acceptors,  that  is,  the  number 
of  acceptors  exceeds  the  shallow  donors.  If  the  number  of 
acceptors  equals  the  number  shallow  donors,  the  middle 
temperature  range  will  be  the  same  as  a  single  donor  plot, 
but  the  low  temperature  range  slope  will  be  equal  to 
average  of  the  activation  energies  of  the  donors.  Only  if 
the  shallowest  donor  is  undercompensated  will  its  activa¬ 
tion  energy  be  given  by  the  slope.  In  all  cases  only  the 
deeper  donors,  which  are  activated  at  higher  temperatures, 
are  unambiguous.  The  concentration  of  the  donor  species 
are  found  directly  from  the  inflection  points  on  the  curve 
(19:13-27). 

Computer  Generated  Curve  Fitting 

Generating  a  fitting  function  to  the  carrier  concen¬ 
tration  data  is  an  iterative  process.  The  parameters, 
acceptor  concentration,  donor  concentrations,  donor  energy 
levels,  and  degeneracy,  are  specified  by  the  user.  The 
computer  generates  the  fitting  function  and  compares  the 
value  of  the  function  to  each  data  point.  It  then  adjusts 
the  parameters,  generates  another  fitting  function,  and 
compares  the  new  values.  This  process  continues  until  the 


user  is  satisfied  that  the  best  fit  has  been  obtained.  The 
computer  reports  the  final  values  of  the  parameters: 
degeneracy,  the  activation  energy,  the  concentration  of 
donors  at  each  energy  level,  and  the  acceptor 
concentration.  The  parameters  may  be  specified  as  fixed  or 
variable.  The  degeneracy  was  fixed  at  two  for  each  donor 
level  in  these  samples.  The  computer  draws  the  fitting 
function  on  the  carrier  concentration  graph  and  gives  the 
percent  error  between  the  function  and  the  data  at  each 
point. 

The  quality  of  the  fit  is  determined  by  the  error 
between  the  data  and  the  fitting  function.  It  is  expected 
that  if  there  are  n  free  parameters,  the  function  should  go 
through  n  points  and  miss  the  rest  by  one  standard 
deviation  (19:28-29).  A  number  value  for  the  quality  of 
the  fit  is  called  Chi-squared  (X2): 

N 

X2  =  £lYi  -  f(Xi)J2/(  <7i)2  (20) 

where 

N  =  number  of  data  points 
Yi  =  i-th  data  point 

ai  =  standard  deviation  at  the  i^h  point 
f(x^)  =  value  of  the  fitting  function  at  x^  (19:30) 

The  standard  deviation  is  assumed  .to  be  0.03Yi  (19:30). 

An  estimate  of  the  minimum  value  of  Chi-squared  depends 
on  the  number  of  data  points,  N,  and  the  number  of  free 
parameters,  n,  so  that  Chi-squared  is  approximately  equal 
to  N  minus  n.  The  quality  of  the  fit  is  then  determined  by 


the  reduced  Chi -squared  (19:30): 

X2V  =  X2/(N-n)  (21  ; 

The  reduced  Chi-squared  is  used  so  that  the  quality  of  a 
one  level  fit  can  be  compared  to  a  multi-level  fit  (23).  A 
multi-level  fit  has  more  free  parameters  than  a  one  level 
fit,  making  the  value  of  Chi-squared  even  smaller  for  the 
same  data  points. 

When  curve  fitting,  it  is  important  to  be  aware  that 
the  program  can  settle  on  unrealistic  values.  This  may  be 
due  to  several  factors  such  as:  systematic  errors  in  the 
experiment,  not  taking  enough  data  points,  or  an  r-f actor 
problem.  If  the  parameter  values  are  unrealistic,  the  user 
can  try  again  with  different  starting  parameters. 


IV.  Experimental  Procedure 

A  boule  of  two  inch  diameter  nitrogen  doped  <100>  FZ 
silicon  manufactured  by  the  Shin-etsu  Company  was  given  to 
the  United  States  Air  Force  by  T.  Abe.  Photoluminescence 
studies  were  performed  by  David  Brown  of  University  of 
Dayton  Research  Institute  for  the  Materials  Laboratory 
( AFWAL ) ,  while  the  Hall  effect  analysis  was  undertaken  in 
this  thesis.  Samples  cut  from  bulk  material  were  annealed 
at  different  temperatures  for  one  hour  then  subjected  to 
the  Hall  experiment  to  determine  the  effect  of  annealing. 
Further,  to  show  the  effect  of  Si:N  as  a  substrate  material 
with  increased  n-type  doping,  a  section  of  the  boule  was 
neutron  transmutation  doped.  The  Hall  effect  experiment 
was  performed  on  samples  from  the  NTD  section  of  Si:N. 

Sample  Geometry  and  Configuration 

In  order  to  accurately  compute  the  Hall  coefficient  and 
resistivity,  the  sample  must  be  homogeneous,  of  uniform 
thickness,  and  contain  no  holes.  All  the  samples  in  these 
experiments  were  square  in  shape  with  contacts  at  each 
corner.  This  configuration  was  suggested  by  L.  J.  van  der 
Pauw  (41).  Van  der  Pauw  showed  that  resistivity  can  be 
accurately  measured  for  an  arbitrarily  shaped  sample  with 
four  point  contacts  on  the  periphery  if  the  above 
conditions  hold  and  that  the  contacts  are  small  compared  to 
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t  h e  ( j'j  t  i  p  1 1  ( ■  t  y  (41:2). 


In  this  conf icjuration  each  contact  plays  a  role  as 
voltage  probe  or  a  current  connection  as  illustrated  in 
Fig.  16.  Measurements  are  required  for  only  three  ol  the 
six  configurations.  Resistivity  measurements  are  taken 
over  two  of  the  four  resistivity  configurations,  1  and  2  oi 
3  and  4.  In  each  configuration  two  voltage  measurements 
are  taken  with  the  current  reversed  for  the  second 
measurement.  Thermoelectric  effects  result  from  the 
dissimilar  metals  used  in  lead  materials  and  are 
independent  of  current  direction.  Averaging  the  resistance 
measurements  eliminates  the  thermoelectric  effects 
( 7: 1 7  ;  27  :  29-31  )  .  The  resistivity  equation  given  by  van  dor 
Pauw  is  (42:221): 

pr  =  lit  (R-j  +  Ra)  f(R-|/R2)  (ohm-cm)  (22) 

"( In 2 )  2 - 

where 

t  =  thickness  (cm) 

Rl  =  V-|  /I  ohms 
R2  =  V2/I  ohms 

The  function,  f(R-|/R2>,  is  a  dimensionless  function  of 
the  ratio  of  the  resistances  of  the  two  separate 
configurations.  If  the  difference  between  R-|  and  R2  is 
small  then  the  function  can  be  approximated  (41:6): 
x  =  R1/R2 

f ( x )  =  1  -  ( x  —  1  2ln(2) 

(x  +  1  H  2 

-  j[  ln(2)  ]  2/4- [ In ( 2 )  )3/12l  (23) 

(x+1  r  1  J 

For  u  symmetrical  sample,  such  as  a  square,  R-| /R2 
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Figure  16.  Van  der  Pauw  Configurations  for 
Hall  Measurements  (19:8) 
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should  be  unity.  A  slightly  irregular  sample  will  have 
some  other  value  for  the  resistance  ratio.  This  ratio 
should  be  nearly  constant  at  all  temperatures.  If  the 
ratio  has  a  temperature  dependence,  it  implies  that  the 
material  is  inhomogeneous  (19:11). 

The  Hall  coefficient,  Rg,  is  obtained  from  either 
configurations  5  or  6.  Measurements  are  taken  over  all 
combinations  of  current  and  magnetic  field  direction.  Four 
voltage  measurements  are  taken  in  all.  Averaging  over 
current  direction  eliminates  thermoelectric  effects,  while 
averaging  over  magnetic  field  direction  eliminates  any 
magneto-resistance  component  due  to  non-synunetric  voltage 
contacts  (7:17).  rh  is  computed  from: 

Z  =  1/2{[R5(  +  ,0+R5<-,*)]/2 

-  [r5(M>+r5<-'H>i/2} 

Rh  =  Zt/B  (cm3/coul)  (24) 

The  plus  or  minus  signs  given  in  parenthesis  after  R5 
refer  to  the  current  direction,  while  the  arrows  refer  to 
the  magnetic  field  direction. 

The  Hall  Effect  Systems 

The  equipment  used  in  these  experiments  was  developed 
by  P.  Hemenger  (9).  This  system  overcomes  the  noise 
problem  and  long  time  constants  encountered  in  measuring 
high  resistivity  samples.  By  using  high-input  impedance, 
unity  gain  amplifiers  between  the  voltage  probe  and  the 
voltmeter  (Hewlett-Packard  3455A  DVM  or  equivalent),  the 


impedance  is  reduced  to  a  level  below  the  input  impedance 
of  the  voltmeter.  In  addition,  by  using  feedback  from  the 
unity  gain  amplifiers  to  drive  the  inner  shield  of  the 
triax  cables  from  the  sample,  the  cable  capacitance  is 
reduced  to  near  zero.  This  cancels  most  of  the 
resistance-capacitance  (RC)  time  constant  which  reduces 
noise,  leakage  current,  and  the  settling  time.  The  system 
uses  Keithley  Model  61 OC  Electrometers  as  the  unity  gain 
amplifiers  which  are  labeled,  A,  in  Figure  17.  A 
laboratory  built  switching  box  uses  reed  switches  to  step 
through  all  the  configurations.  This  system  is  capable  of 
measuring  resistances  up  to  1012  ohms. 

The  Hall  effect  system  is  schematically  shown  in  Figure 
18.  Data  was  collected  on  two  similar  systems.  The 
function  of  each  block  is  identical  for  both  systems.  In 
each  system  the  current  is  monitored  by  an  ammeter  capable 
of  measuring  current  in  the  order  of  10~^  amps  (19:5). 

The  system  referred  to  as  Hall  I  can  be  operated  manually 
or  by  computer  (ADAC  Corp.  LSI  System  100).  Hall  II  can 
only  be  operated  manually.  The  other  differences  between 
the  systems  are  that  Hall  II  uses  a  constant  current  source 
(Keithley  220  Programmable  Current  Source)  where  Hall  I 
uses  a  DC  voltage  power  supply  (Hewlett-Packard  6299A)  to 
provide  sample  current  and  that  the  gaussmeters  are  made  by 
different  manufacturers. 

The  computer  that  is  used  to  automatically  operate  Hall 
I  is  not  shown  in  Figure  18.  The  computer  controls  the 


Hall  System 


Sample  Temperature  Magnet 


temperature,  magnet,  the  configuration  of  the  sample  by 
setting  the  reed  switches  on  the  Van  der  Pauw  Box,  and 
collecting  the  data  by  reading  the  meters.  When  operated 
automatically,  Hall  I  collects  data  over  all  resistivity 
and  Hall  configurations,  taking  16  measurements  in  all. 

The  temperature  of  the  sample  is  determined  by 
measuring  the  voltage  drop  across  a  reverse  biased  silicon 
diode  thermometer  (Lake  Shore  Cryogenics  Model  DT500).  The 
voltage  drop  across  the  diode  is  read  by  a  Hewlett-Packard 
3455A  Digital  Voltmeter.  The  temperature  is  determined  by 
referencing  the  voltage/temperature  calibration  data 
provided  by  the  manufacturer. 

The  sample  probe  is  mounted  in  a  Janis  Super  Veri-Temp 
helium  vapor  dewar.  The  dewar  is  filled  with  liquid 
helium.  The  helium  flows  through  a  needle  valve  at  the 
bottom  of  the  dewar  then  vaporizes.  This  allows  cold 
helium  vapors  to  flow  around  the  sample  holder.  The  sample 
holder  (Fig.  19)  is  a  block  of  copper  which  acts  as  a  heat 
sink  and  is  covered  by  a  thin  cylinder  of  copper  as  a 
radiation  shield.  The  sample  is  mounted  on  a  beryllium 
oxide  (BeO)  pad  and  held  in  place  by  a  non-magnetic  spring 
with  a  thin  3heet  of  teflon  in  between  for  electrical 
insulation.  BeO  electrically  insulates  the  sample  from  the 
copper  sample  holder  while  keeping  in  thermal  contact  with 
the  silicon  diode  thermometer.  BeO  is  an  unusual  material 
that  has  high  thermal  conductivity  but  very  low  electrical 
conductivity.  The  silicon  diode  thermometer  is  directly 
underneath  and  in  contact  with  the  BeO  pad. 


Figure  19.  Schematic  Diagram  of  Sample  Probe 
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To  raise  the  sample  temperature  above  liquid  helium 


temperature,  two  heater  coils  are  positioned  at  each  end  of 
the  sample  holder.  The  heater  coils  are  wound  in  opposite 
directions  to  reduce  maqnet.ic  fields  (7:29).  When  takinq 
data  manually  on  either  system,  sample  temperature  is 
maintained  by  an  Artronix  Instrumentation  Model  5301 
temperature  controller.  The  controller  senses  temperature' 
through  one  of  the  two  thermometers  near  the  bottom  of  the 
sample  holder.  The  carbon-glass  thermometer  used  for 
temperatures  between  4  to  30  K  while  the  platinum 
thermometer  is  used  for  temperatures  30  K  and  above.  These 
thermometers  unlike  the  Si  diode  thermometer  are 
insensitive  to  magnetic  fields  in  these  temperature 
regions.  The  temperature  controller  drives  the  heater 
coils.  When  automatically  taking  data.  Hall  I  uses  a 
laboratory  built  temperature  controller  that  is  under 
computer  control.  The  automatic  temperature  controller 
uses  the  same  controlling  thermometers  in  the  sample 
holder. 

The  magnetic  field  is  supplied  by  a  Varian 
electromagnet  which  monitored  by  a  Walker  Magnametric  Model 
MG-3D  gaussmeter  (Hall  I)  or  Bell  Model  620  gaussmeter 
(Hall  II).  The  sample  probe  is  mounted  in  the  dewar  so 
that  the  magnetic  field  is  normal  to  the  sample  surface. 

Sample  Preparation 

After  the  samples  are  cut,  they  must  be  prepared  for 
annealing.  Contacts  must  be  put  on  the  corners  and  then 


mounted  on  the  sample  holder.  Each  step  requires  care  to 
prevent  contamination. 

Prior  to  annealing  the  samples  are  rinsed  successively 
in  trichloroethylene,  acetone,  then  methanol.  Next  they 
are  placed  in  a  solution  of  ammonium  hydroxide,  hydrogen 
peroxide,  and  water  (1:1:5)  at  60  °C  for  15  minutes. 

After  this  they  are  placed  in  a  solution  of  hydrochloric 
acid,  hydrogen  peroxide,  and  water  (1:1:6)  at  60  °C  for 
15  minutes.  The  samples  are  rinsed  in  de-ionized  water  and 
allowed  to  soak  for  30  minutes  or  more. 

The  boat  used  for  annealing  the  samples  is  high  purity 
polycrystalline  silicon.  Wafers  of  the  same  material  are 
placed  between  the  samples  to  prevent  possible 
contamination  during  annealing.  Before  annealing  the  boat 
and  the  wafers  are  cleaned  using  a  similar  procedure. 

During  annealing,  argon  gas  is  allowed  to  flow  through  the 
furnace. 

The  first  two  samples  were  annealed  at  900  °C  using  a 
quartz  furnace  tube.  Unfortunately,  the  tube  broke  during 
cooling.  Fiberglass  was  used  to  pack  the  furnace  tube  to 
keep  it  from  losing  too  much  heat.  The  fiberglass  may  have 
contracted  faster  than  the  quartz  tube  causing  it  to 
fracture  (23).  After  this  a  polycrystalline  silicon 
furnace  tube  was  used,  but  the  temperature  could  not  be  set 
above  830  °C. 

To  make  ohmic  contacts  on  the  corners  of  the  samples,  a 
pulsed  Nd: YAG  laser  is  used  to  melt  a  1  mm  spot  while 


driving  in  dopant  material.  The  samples  were  immersed  in  a 
solution  of  spin-on  phosphorus  dopant.  The  laser  pulse 
melts  a  small  localized  area  which  incorporates  phosphorus 
into  the  lattice  as  it  recrystallizes.  The  result  is  a 
degenerately  doped  area  suitable  for  an  ohmic  contact.  To 
remove  the  phosphorus  spin-on  from  the  surface  of  the 
samples,  they  were  placed  in  a  solution  of  ammonium 
hydroxide,  hydrogen  peroxide,  and  water  (1:1:5)  at  60  °C 
for  15  minutes,  then  rinsed  in  distilled  water. 

After  cleaning,  copper  leads  are  soldered  on  using 
indium  solder.  The  samples  are  mounted  to  the  holder  by 
soldering  the  leads  to  the  center  wire  of  the  shielded 
cables  on  the  sample  probe  with  indium  solder.  The  probe 
is  ready  for  mounting  in  the  Hall  equipment. 

As  a  check  before  each  data  run,  the  sample  is  tested 
for  ohmic  contacts.  A  current  is  applied  to  the  sample  and 
the  resulting  voltages  are  read  for  each  of  the  resistivity 
configurations  1  through  4  in  both  current  directions.  The 
voltage  at  each  configuration  should  be  the  same  for  both 
current  directions.  If  not,  one  or  more  contacts  are  not 
ohmic.  Increasing  the  current  should  result  in  a 
proportional  increase  in  the  voltage.  This  test  will  also 
reveal  broken  or  shorted  wires. 


V.  Permits 


The  Hall  effect  measurements  were  performed  on  three 
different  types  of  silicon:  phosphorus  doped  CZ  silicon, 
nitrogen  doped  FZ  silicon,  and  neutron  transmutated  (NTD) 
nitrogen  doped  FZ  silicon.  All  of  the  nitrogen  doped 
samples,  including  the  NTH  samples,  came  from  the  same 
boule.  The  CZ  was  used  as  an  example  for  comparison  with 
the  FZ  Si:N  samples,  while  a  high  purity  FZ  NT1)  sample 
(NT2-07-5N)  that  had  been  analyzed  by  W.  C.  Mitchel  was 
used  to  compare  results  from  the  NTD  Si:N  samples  (23). 

The  Si:N  samples  were  annealed  for  one  hour  at  different 
temperatures  to  show  the  effects  of  annealing.  Table  II 
lists  the  sample  identification  number,  silicon  type,  dop¬ 
ing,  and  annealing  temperature  as  a  summary  of  the  exper¬ 
iments  performed.  The  magnetic  field  strength  was  0.50 
Tesla  for  all  experiments  except  FZ0001  which  was  0.10 
Tesla.  Due  to  surface  contamination,  sample  0274-0802  and 
0274-802-1426  were  recleaned  in  the  ammonium  hydroxide, 
hydrogen  peroxide,  and  water  solution  and  the  experiments 
repeated. 

Electrical  Data  at  300  K 

Since  most  electrical  devices  are  designed  to  operate 
at  room  temperature,  Table  III  gives  the  Hall  coefficient, 
Rjj ;  resistivity,  pr;  and  Hall  mobility,  uh»  at 
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V 


Table  II 


List  of  Samples  on  Which  Experiments  Were  Performed 


Name 

in# 

Type 

Doping 

Element 

Annealing 

Temperature 

CZ-AN-Q01 

01 24-1 399-131 2 

CZ 

P 

900  °C 

FZ0002 

0274A-1 581 -1 395 

FZ 

N 

800  °C 

( NTD ) 

P 

FZ0003 

0274A-1 581 -1 396 

FZ 

N 

800  °C 

(NTD) 

P 

FZ0001 

0274-0802 

FZ 

N 

900  °C 

FZ0001 -A 

0274-0802 

FZ 

N 

900  o c 

FZ0004 

0274-802-1 426 

FZ 

N 

800  °C 

FZ0004-A 

0274-802-1426 

FZ 

N 

800  °C 

FZ0005 

0274-802-1422 

FZ 

N 

830  °C 

FZ0006 

0274-802-1427 

FZ 

N 

unannealed 

SI-N800C 

0274-802-1 427 

FZ 

N 

800  oc 

Table  III 


Electrical  Parameters  at  300  K 


Sample 
&  Temp. 

l1/ 

cm-3/ coul 

Pr 

ohm -cm 

Carri er 
Cone.  cm-3 

9  UH 

cnr/V-sec 

Annealed  Si 

:N  Samples 

FZ0006 
unannca led 

1.62  X  1 O6 

1  1  80 

3.85  X  1012 

1  380 

SI-N800C 

800  °C 

2.10  X  105 

302 

2.97  X  1013 

696 

FZ0004-A 

800  °C 

8.60  X  105 

652 

7.26  X  1012 

1  320 

FZ0005 

830  °C 

4.76  X  105 

801 

1.31  X  1013 

593 

FZ0001 -A 

900  °C 

6.49  X  105 

600 

9.62  X  1012 

1080 

Phosphorus 

Doped  CZ  Sample 

CZ-AN-001 
900  °C 

2.20  X  103 

1.62 

2.84  X  1015 

1  360 

NTD  Si:N  Samples 

6.88  X  104 

35.6 

9.07  X  1013 

4.66  X  105 

34.1 

1.34  X  1014 

FZ0002 
800  °C 

FZ0003 
800  °C 


High  Purity  NTD  FZ  Sample 


1  930 

1  370 


NT2-07-SN 
800  °C 


1.95  X  1 05 


1  281 


1.34  X  1014 


1  530 


PV'J  .«  J 


300  Kelvin.  The  results  for  the  phosphorus  doped  sample, 
CZ-AN-001 ,  are  within  5.3%  of  the  expected  results  for 
resistivity  and  4.6%  for  mobility  when  compared  to  the  data 
in  Figure  20  and  21  (37:29,32).  It  should  be  recalled  that 
although  drift  mobility  is  equal  to  conductivity  mobility 
where  one  conduction  species  (either  electrons  or  holes) 
dominates,  the  Hall  mobility  is  equal  to  the  conductivity 
mobility  only  when  the  r-factor,  r,  is  equal  to  one.  Thus, 
conclusions  about  conductivity  mobility  derived  from  the 
Hall  mobility  data  should  not  be  made  too  quickly. 

Assuming  for  the  moment  that  in  these  experiments  r  = 

1 ,  Table  IV  gives  a  comparison  of  the  experimental  values 
for  carrier  concentration  and  mobility  with  values  taken 
from  the  graphs  in  Figures  20  and  21 .  Resistivity  data  was 
used  to  determine  the  carrier  concentration  because  it  is 
directly  measured.  Most  of  the  samples  have  carrier 
concentrations  in  the  range  of  lO1^  to  10^  3 
electrons/cm3  which  is  below  the  range  given  in  Fig.  21 . 

At  low  concentrations  impurity  scattering  is  negligible  and 
the  mobility  is  entirely  determined  by  phonon  scattering 
( 21 :308-31 7;23 ) .  From  the  graph  of  Figure  21  this  limit  is 
reached  around  an  impurity  concentration  of  lO1^  cm-3. 

The  limiting  value  is  1500  cm^/v-sec.  However,  the 
mobility  of  the  Si:N  NTD  samples  vary  as  much  as  28.6% 
while  the  high  purity  sample  (NT2-07-SN)  varies  by  8.7%. 

The  Si:N  samples  should  not  be  compared  with  the 
resistivity /concentration  graph  in  Fig.  20,  as  this  graph 
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Figure  20.  Resistivity  versus  Impurity  Concentration 
for  Silicon  at  300  K  (37:32) 


IMPURITY  CONCFNTRATION  (cm  3) 


Figure  21 .  Drift  Mobility  of  Silicon  versus 

Impurity  Concentration  at  300  K  (37:29) 
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Ta  bio  IV 


Comparison  of  Experimental  Data  with  Values  from 

Figures  20  and  21 


Sample 

Pr 

ohm -cm 

Carrier 

Cone.  cm-3 

Mobility 

cm3/V-sec 

Experiment:  CZ-AN-001 

1  .62 

2.84  X  1015 

1  360 

Sze's  data: 

3  X  101s 

1  500 

Deviation : 

5.3% 

9.3% 

Experiment:  FZ0002 

35.6 

9.07  X  1 01 3 

1930 

Sze's  data: 

9  X  1013 

1  500 

Deviation : 

0.7% 

28.6% 

Experiment:  FZ0003 

34.1 

13.4  X  101 3 

1  3  70 

Sze's  data: 

13  X  1013 

1  500 

Deviation : 

3.1  % 

6.7% 

Experiment:  NT2-07-SN 

1  28 

3.21  X  1013 

1  520 

Sze's  data: 

3  X  1 01  3 

1  500 

Deviation : 

7.0% 

1  .3% 

was  generated  from  phosphorus  doped  silicon  data.  The 
generalized  resistivity  curves  for  n  and  p  type  silicon 
published  by  Irvin  gives  data  for  material  doped  in  the 
range  of  10^  to  lO^1  cm-3  (14).  Irvin's  curve  for 
n-type  material  was  compiled  from  various  Four-point  probe 
resistivity  and  Hall  experiments  performed  on  silicon 
material  doped  with  group  V  donors,  but  not  nitrogen 
(14:388-391).  However,  if  the  electron  mobility  is  assumed 
to  be  about  1500  cm^/v-sec,  then  the  data  in  Table  III 
shows  that  the  Si:N  samples  have  much  lower  than  expected 
mobilities.  Resistivity  and  mobility  is  the  highest  lor 
the  unannealed  sample,  FZ0006,  but  its  concentration  is  the 
lowest.  This  is  consistent  with  Fig.  21,  which  shows  that 
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mobility  decreases  with  increasing  concentration.  At  1380 
cm^/v-sec  the  mobility  of  FZ0006  is  nearly  the  same  as 
CZ-AN-001  at  1360  cm^/v-sec,  but  the  concentration 
differs  by  1000.  After  FZ0006  was  annealed,  the  sample 
name  was  changed  to  SI-N800C  in  order  to  keep  the  data 
files  separate.  Table  III  shows  that  resistivity,  carrier 
concentration,  and  mobility  all  decreased.  Although  Table 
III  lists  Si:N  samples  with  increasing  annealing 
temperatures,  no  clear  trend  is  established,  but  it  appears 
that  annealing  does  affect  the  material. 

Comparison  of  Electrical  Data  from  20  to  380  K 

The  Hall  experiments  were  performed  at  temperatures 
ranging  from  20  to  380  K.  The  plots  of  temperature 
dependent  data  give  more  information  as  to  the  nature  of 
the  material.  Comparisons  are  made  by  plotting  various 
samples  on  the  same  graph.  Plots  were  made  of  the  carrier 
concentration,  Hall  coefficient,  resistivity,  and  Hall 
mobility  versus  temperature. 

Figures  22  through  23  plot  the  data  for  the  Si:N 
samples:  FZ0001-A,  FZ0006,  and  SI-N800C.  FZ0006  is  the 

unannealed  Si;N  sample,  while  SI-N800C  is  the  same  sample 
annealed  at  800°C.  FZ0001-A  was  annealed  at  900°C 

Data  below  40  K  could  not  be  obtained  for  FZ0006. 
Resistivity  above  10^  ohm-cm,  resulted  in  excessively 


long  settling  time. 

The  Hall  coefficient,  which  is  a  directly  measured 


c? 


Figure  22.  Hall  Coefficient  of  Three  Si:N  Samples 
versus  Inverse  Temperature 
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Figure  25.  Hall  Mobility  of  Three  Si:N  Samples  versus 
Inverse  Temperature 


parameter,  is  plotted  with  inverse  temperature  (Fig  22). 
After  annealing  at  800°C,  the  unannealed  sample  (FZ0006) 
has  changed  considerably  in  the  high  temperature  region 
(SI-N800C).  It  does  not  track  the  unannealed  FZ0006  or  the 
900°C  sample  (FZ0001-A).  FZ0006  and  FZ001 -A  track  each 

other  well  in  the  low  temperature  regions  but  start  to 
separate  above  106  K. 

The  other  directly  measured  parameter  is  resistivity. 
Resistivity  is  also  plotted  against  inverse  temperature 
(Fig  23).  The  plot  of  SI-N800C  appears  to  be  highly 
irregular  compared  to  the  other  samples  and  FZ0006.  In 
fact,  FZ0001 -A  and  FZ0006  are  typical  resistivity  plots. 
Above  55  K,  in  the  resistivity  plots  of  FZ0001 -A  and 
FZ0006,  resistivity  increases  with  temperature.  At  this 
temperature  all  of  the  donors  are  ionized  and  the  increase 
in  resistivity  is  from  phonon  scattering  (21:277-281, 
308-317).  Below  55  K,  carrier  concentration  decreases  so 
resistivity  increases.  The  concentration  versus  inverse 
temperature  plot.  Figure  24,  shows  that  the  carrier 
concentration  is  saturated  at  55  K. 

Carrier  concentration  which  is  computed  from  Eq.  9 
assumes  an  r-factor  of  one.  Thus,  the  carrier 
concentration  is  inversely  proportional  to  Rjj.  The 
concentration  versus  temperature  curve  in  Fig.  24  shows 
that  the  carrier  concentration  levels  off  at  high 
temperatures  for  the  unannealed  sample,  FZ0006,  but 
annealing  causes  a  change  in  the  material  so  that 
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concentration  starts  to  increase  again.  FZ0001-A  exhibits 
the  same  behavior  in  the  high  temperature  region.  This 
implies  that  deep  levels  are  activated  by  high  temperature 
annealing,  but  it  appears  that  increasing  the  temperature 
does  not  result  in  the  proportional  increase  in  carriers. 

For  the  readers  convenience.  Hall  mobility  is  plotted 
both  with  inverse  temperature  and  temperature  in  Figures  25 
and  26.  Hall  mobility  varies  proportionally  with  and 
inversely  with  pr  (Eq.  11).  For  FZ0001 -A  mobility  peaks 
at  about  7  x  1 04  cm2/V-s ec  at  30  K.  SI-N800C  has  a 
peculiar  decrease  and  increase  below  85  K  and  and  is 
displaced  upwards  at  225  K.  This  is  actually  not  so 
peculiar  in  that  the  resistivity  and  the  Hall  coefficient 
curves  were  not  typical.  What  is  peculiar  is  that  a  sample 
that  produced  typical  data  has  an  anomalous  pr  after 
annealing,  but  another  sample  annealed  at  a  higher 
temperature  behaves  quite  well. 

The  next  set  of  figures  show  the  phosphorus  doped  CZ 
sample  is  plotted  with  an  FZ  sample,  FZ0001 -A.  These  were 
both  annealed  at  the  same  time  at  900  °C.  The  Hall 
coefficient  plot,  Fig.  27,  shows  that  the  Rfj  of  the 
phosphorus  doped  CZ  sample,  CZ-AN-001 ,  decreases  slightly 
as  the  temperature  decreases  in  the  380  to  200  K  region 
then  starts  to  increase  as  the  temperature  decreases. 
However,  Fig.  28  shows  a  well  behaved  resistivity  plot. 

The  carrier  concentration  plot.  Fig.  29,  shows  the  carrier 
concentration  of  CZ-AN-001  decreasing  at  above  200  K. 
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Figure  27.  Hall  Coefficient  of  CZ  and  Si:N  Samples 
Annealed  at  900  °C 
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Figure  28.  Resistivity  of  CZ  and  Si:N  Samples 
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Figure  30.  Hall  Mobility  of  CZ  and  Si:N  Sample 
Annealed  at  900  °C  versus  Inverse  Temperature 


The  dip  in  in  at  the  200  to  300  K  region  is  due  to 
temperature  dependence  in  the  scattering  mechanism. 
Scattering  effects  are  included  in  the  r-factor  which  was 
assumed  to  be  unity  at  all  temperatures.  Although  very 
difficult  to  compute,  the  r-factor  is  known  to  have  some 
temperature  dependence  (19:66,67).  This  also  explains  the 
decrease  in  carrier  concentration  above  200  K.  The  carrier 
concentration  is  calculated  from  Eq.  9,  which  is  restated 
here : 


n  =  r 

( cm -3 ) 

(9) 

eRH 

becomes 

n  =  r ( T ) n0 

( cnr^  ) 

(25) 

where  nQ  is  the  experimental  value  of  the  concentration 
and  r(T)  is  the  temperature  dependent  r-factor.  A  careful 
re-examination  of  FZ0006  in  Fig.  24  shows  a  slight  decrease 
in  concentration  at  high  temperatures.  However,  the 
activation  of  deeper  energy  levels  in  Si:N  by  annealing 
seems  to  overcome  the  pulling  down  of  the  concentration  by 
the  r-factor. 

Hall  mobility  is  computed  using  Eq.  11  and  the  r-factor 
appears  in  the  Hall  coefficient,  but  the  effect  of  the 
r-factor  is  not  clear  in  the  Hall  mobility  plots  (Fig.  30 
and  31)  of  CZ-AN-001.  The  decrease  in  may  be  offset 
by  the  resistivity.  It  is  interesting  to  note  that 
CZ-AN-001  has  higher  mobility  than  FZ0001-A,  but  only  above 
200  K.  FZ0001 -A  has  higher  Hall  mobility  below  this 
temperature.  However,  because  of  the  r-factor  problem,  it 
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is  unwise  to  say  that  the  conduction  mobility  is  exactly 
equal  to  the  Hall  mobility. 

Combining  the  results  of  the  NTD  samples  into  one  set 
of  graphs  (Figures  32  to  36)  show  the  Si:N  samples  differ 
from  the  high  purity  sample  (NT2-07-5N)  mainly  in  the 
amount  of  phosphorus  doping.  All  three  track  each  other 
well  in  the  Hall  coefficient  (Fig.  32)  and  resistivity 
(Fig.  33)  graphs.  The  figures  show  that  FZ0002  is  not  as 
well  behaved  as  the  other  samples.  FZ0002  was  performed 
manually  on  Hall  II  while  FZ0003  and  NT2-07-5N  were 
performed  automatically  on  Hall  I.  Thus,  the  slightly 
erratic  behavior  of  FZ0002  may  be  due  to  increased 
experimental  uncertainty.  Again  there  is  a  slight  increase 
in  Rh  at  high  temperatures  in  all  three  plots  indicating 
an  r-factor  problem.  The  concentration  graph.  Figure  34, 
also  has  the  same  decrease  in  concentration  at  high 
temperatures  as  CZ-AN-001  for  all  three  plots.  However, 
the  Hall  mobility  graphs  (Figures  35  and  36)  show  that 
FZ0002  and  FZ0003  have  the  same  type  behavior  at  low 
temperatures  as  seen  in  SI-N800C  (Fig.  25  and  26). 

The  carrier  concentration  graph  (Fig.  34)  shows  that  at 
very  low  temperatures,  less  than  28  K,  the  FZ0002  and 
FZ0003  plots  appear  to  bend  away  from  the  NT2-07-5N  plot  on 
the  concentration  graph.  This  bending  outward  at  low 
temperatures  also  appears  in  Fig.  24  in  the  plot  of 
FZ0001-A.  This  may  indicate  an  energy  level  shallower  than 
phosphorus  in  the  Si:N  material.  However,  there  are  other 
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Figure  36.  Hall  Mobility  of  si:N  NTD  and  High  Purity  FZ 
NTD  Samples  versus  Temperature 


parameters  to  be  considered  before  any  statement  can  be 
made. 

Temperature  Dependence  of  Ri /R? 

In  Chapter  IV  it  was  shown  that  resistivity  as  given  by 
van  der  Pauw  (Eq.  22)  has  an  f(R-|/R2)  term  which  should 
be  constant  over  all  temperatures,  if  not,  inhomogeneity  is 
implied.  This  parameter  is  calculated  for  every  data  point 
in  the  experiment.  It  is  not  unusual  to  see  some 
temperature  dependence  at  low  temperatures. 

This  was  the  case  for  CZ-AN-001 ,  the  phosphorus  doped 
CZ  sample.  Below  30  K  there  is  some  temperature  dependence 
in  R1/R2  as  seen  in  Fig.  37.  Figures  38  through  40 
plot  the  CZ-AN-001  data  from  30  to  380  K;  however,  the 
experimental  data  started  at  20  K.  Figures  38  and  39  show 
a  decrease  in  the  Hall  coefficient  and  an  increase  in 
concentration  below  30  K.  This  is  believed  to  be  due  to 
surface  conduction  across  the  sample,  because  the 
resistivity  does  not  go  above  400  ohm-cm,  as  seen  in  Fig. 
39.  This  phenomenon  disappears  at  30  K  as  conduction 
through  the  sample  takes  over.  The  simplest  way  to  deal 
with  this  is  to  eliminate  the  data  below  28  K.  In  any 
event,  data  in  this  region  should  not  be  relied  on.  This 
is  why  Figures  27  through  31  were  plotted  without  the  low 
temperature  points. 

If  temperature  dependence  is  seen  over  the  entire 
range,  then  the  sample  is  inhomogenous.  Thus,  a  necessary 
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Figure  37.  R-j  /R2  Ratio  of  CZ-AN-001 
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Figure  39.  Resistivity  of  CZ-AN-001  from  20  -  380  K 
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condition  of  the  experiment  has  failed.  Not  only  is  the 
resistivity  data  questionable,  but  the  Hall  coefficient 
data  is  also  unreliable  due  to  the  unknown  effect  of 
inhomogeneity  on  the  scattering  mechanism  (r-factor). 
Unfortunately,  this  was  the  case  for  FZ0004-A  and  FZ0005 
seen  in  Figures  41  and  42.  Sample  0274-802-1426 
(experiment  FZ0004  and  FZ0004-A)  was  annealed  at  800  °C 
and  sample  0274-802-1422  (experiment  FZ0005)  was  annealed 
at  830  °C.  The  data  for  the  first  experiment  on 
0274-802-1426  (FZ0004)  showed  a  temperature  dependent 
R-| /R2  and  anomalous  Rjj,  pr,  u^,  and  concentration 
plots.  It  was  thought  that  surface  contamination  was 
responsible  so  the  sample  was  cleaned  in  a  solution  of 
ammonium  hydroxide,  hydrogen  peroxide,  and  water  and  new 
leads  soldered  on.  The  results  of  the  second  experiment 
were  no  better.  FZ0005  also  had  a  temperature  dependent 
R-|/R2.  The  graphs  of  the  data  for  FZ0004-A  and  FZ0005 
are  given  in  Appendix  B. 

Of  the  Si:N  samples,  FZ0001-A  and  FZ0006  have 
acceptable  R-|  /R2  ratios,  but  SI-N800C  does  not.  FZ0006 
and  SI-N800C  are  the  same  sample.  Before  it  was  annealed 
FZ0006  had  an  acceptable  R1/R2  ratio  (Figure  43). 

Figure  44  shows  that  after  annealing  the  sample  (renamed 
SI-N800C)  is  temperature  dependent  over  its  entire  range 
and  that  the  difference  between  the  highest  and  lowest 
points  is  0.36.  This  indicates  that  annealing  produced  a 
change  in  the  sample  that  caused  it  to  become 
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Figure  41.  H1/R2  Ratio  of  Si:N  Sample  Annealed 
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Figure  43.  Ri/R?  Ratio  of  Unannealed  Si:N  Sample 
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inhomogeneous.  FZ0001 -A  has  a  pronounced  temperature 
dependence  in  Ri /R2  at  low  temperatures  as  seen  in  Fig. 

45.  If  the  points  below  27  K  are  eliminated,  then  the 
difference  between  the  highest  and  lowest  point  is  0.16. 

Some  temperature  dependence  in  R1/R2  *s  n°l 
surprising  in  the  NTD  samples.  The  NTD  process  damages  the 
lattice;  however,  almost  all  of  the  damage  is  repaired  by 
annealing  (23).  The  residual  damage  results  in  slightly 
inhomogeneous  material.  Figure  46  shows  the  R1/R2  plot 
for  NT2-07-5N.  The  ratio  is  temperature  dependent  but  has 
a  variation  of  0.136.  The  R;/R2  ratio  in  FZ0002 
appears  to  have  a  strong  temperature  dependence  in 
R1/R2  below  45  K  and  and  above  45  K  the  R1/R2  plot 
seems  to  become  a  straight  line  (Figure  47).  However,  the 
variation  of  R1/R2  at  temperatures  from  45  to  380  K 
have  a  variation  of  0.4.  This  experiment  was  performed 
manually  and  may  have  more  experimental  uncertainty.  The 
data  points  display  a  random  pattern  in  this  region  rather 
than  a  smooth  upwards  or  downward  trend.  Thus,  this 
R-| /R2  is  marginally  acceptable  above  40  K.  The  last 
sample,  FZ0003,  also  shows  a  strong  temperature  dependence 
at  low  temperatures  and  then  levels  out  (Fig.  48).  The 
total  variation  in  R-j /R2  is  0.80  over  the  entire 
temperature  range,  but  above  40  K  the  variation  is  0.12. 
This  means  that  the  R1/R2  in  FZ0003  is  acceptable  at 
temperatures  above  40  K. 
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Figure  45.  Ri/Ro  Ratio  of  Si:N  Sample  Annealed 
at  900  °C 
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Figure  48.  Ri /R2  Ratio  of  Si:N  NTD  Sample  (FZ0003) 


Data 


Curve  Fitting  of  Experimental 

Attempts  to  use  the  curve  fitting  computer  program  were 
hindered  by  the  r-factor  problem  and  the  temperature 
dependent  R1/R2  ratio  at  low  temperatures  in  the 
samples.  Because  the  data  is  unreliable  where  R-| /R2  is 
temperature  dependent,  these  data  points  are  eliminated. 

The  r-factor  problem  effects  every  data  point,  but  has  the 
most  obvious  effect  at  high  temperatures.  Eliminating  the 
points  where  the  carrier  concentration  is  decreasing  with 
temperature  helps  generate  a  better  fit.  Nevertheless, 
there  will  still  be  an  error  in  the  activation  energy  and 
the  concentration  of  tne  donors  (19:36-38).  Due  to  the 
r-factor  problem,  if  the  goodness  of  fit  parameter. 
Chi-squared  per  degree  of  freedom,  is  less  than  three  it 
can  be  considered  as  the  best  fit  (23). 

The  phosphorus  doped  CZ  sample,  CZ-AN-001 ,  had  a 
temperature  dependent  R1/R2  at  low  temperatures  and  an 
apparent  r-factor  problem  at  high  temperatures.  Elimin¬ 
ating  points  above  200  K  and  below  30  K  still  left  enough 
data  points  to  generate  a  good  curve  fitting  function.  The 
results  as  seen  in  Fig.  15  give  a  Chi-squared  value  of 
0.39,  which  is  an  excellent  fit.  The  donor  was  identified 
from  the  activation  energy  as  phosphorus  and  the  donor 
concentration  was  given  as  3.19  x  10^5  cm~3.  Because 
of  the  r-factor  the  phosphorus  concentration  should  be 
given  as  about  3  x  10^  cm-3. 

The  NTD  samples  also  exhibited  an  obvious  r-factor 


93 


problem  at  high  temperatures.  Curve  fitting  the  high 
purity  NTD  FZ  sample,  NT2-07-5N,  was  accomplished  by 
eliminating  points  above  140  K.  The  best  fit  had  a 
Chi-squared  of  4.05  as  seen  in  Figure  49.  Table  5  gives 
the  results  for  the  NTD  samples.  For  NT2-07-5N  the  number 
of  donors  is  3.68  x  10^3  cm“3  and  the  activation  energy 
is  0.0456  eV.  This  means  that  the  sample  was  doped  with 
phosphorus  at  a  concentration  of  about  4  x  10^3  cm~3. 

The  curve  fitting  program  cannot  give  any  information  about 
the  number  of  acceptors  other  than  the  concentration.  It 
is  likely  that  the  acceptors  are  boron  atoms.  The 
segregation  coefficient  (ratio  of  the  concentration  in 
solid  silicon  to  the  concentration  in  melted  silicon)  of 
boron  is  close  to  unity  so  that  it  is  virtually  impossible 
to  zone  refine  a  boule  of  silicon  to  have  a  concentration 
of  boron  less  than  10^  cm~3  (23). 

Due  to  the  r-factor  problem,  points  above  120  K  were 
deleted  from  the  Si:N  NTD  samples.  The  R1/R2  ratio  in 
these  samples  was  very  temperature  dependent  at  low 
temperatures,  so  points  below  32  K  were  also  eliminated. 

The  curve  fitting  program  resulted  in  Chi-squared  values  of 
1.38  for  FZ0002  and  3.67  for  FZ0003.  Figures  50  and  51 
give  the  resulting  curves.  Again  the  activation  energy  is 
affected  by  the  r-factor,  but  it  is  close  enough  to  be 
identified  as  phosphorus.  Clearly  the  Si:N  NTD  samples  are 
dominated  by  the  phosphorus  and  behave  similarly  to  the 
high  purity  FZ  NTD  sample. 
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Figure  50.  Least  Squares  Fit  of  Carrier  Concentration 
of  SisN  NTD  Sample  (FZ0002) 


Table  V 

Results  of  Curve  Fitting  NTD  Samples 


Sample 

Acceptor 
Cone,  cm-2 

Donor 

Cone.  cm-3 

Activation 
Energy  eV 

X2 

NT2-07-5N 

1.69  x  1012 

3.68  x  1013 

0.0456 

4.05 

FZ0002 

3.85  x  1012 

.-A 

• 

O 

O' 

X 

o 

0.0485 

1 .38 

FZ0003 

4.64  x  1012 

1  .51  x  1  014 

0.0481 

3.67 

The  concentration  data  was  much  more  difficult  to 
computer  fit  for  the  rest  of  the  Si:N  samples.  It  was 
impossible  to  try  to  fit  FZ0004-A,  FZ0005,  and  SI-N800C 
because  of  the  anomalous  carrier  concentration  plots  and 
more  importantly  because  of  the  temperature  dependence  in 
R-|/R2  over  all  temperatures.  This  left  the  unannealed 
Si:N  sample,  FZ0006,  and  the  900  °C  annealed  Si:N  sample, 
FZ0001-A,  for  curve  fitting. 

A  one  level  fit  was  performed  on  FZ0006,  resulting  in  a 
solution  of  4.58  x  10^  cm-3  acceptors,  4.02  x  1012 
cm-3  donors,  and  an  activation  energy  of  0.053  eV. 
Chi-squared  was  3.48.  Next,  a  two  level  fit  was  attempted 
in  order  to  generate  a  better  fit  (Figure  52).  The  deeper 
energy  level  turned  out  to  be  0.50  eV  at  2.13  x  1012 
cm-3;  however,  this  level  is  probably  indicates  the  onset 
of  intrinsic  conduction.  At  high  temperatures  there  is 
enough  thermal  energy  to  start  ionizing  silicon  atoms 
( 27: 1 06-1 09 ) .  Thus  a  sharp  rise  in  concentrat ion  is  seen 
at  high  temperatures.  The  lower  energy  level  was  0.052  eV 
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Figure  52.  Two  Level  Least  Squares  Fit  of  Carrier 

Concentration  of  Unannealed  Si:N  Sample 
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at  3.96  x  10  cm  and  the  acceptor  concentration  was 
3.41  x  10^0  cm  ~  3 . 

The  only  sample  with  an  acceptable  R1/R2  ratio 
after  annealing  was  FZ0001-A.  Points  below  28  °K  were 
eliminated  as  R1/R2  was  temperature  dependent  in  this 
region.  The  best  one  level  fit  on  data  from  28  to  380  K 
gave  a  Chi-squared  of  73.48  (Fig.  53).  Most  of  the  error 
appears  to  be  in  the  high  temperature  region.  In  this  case 
it  can  be  said  that  deep  levels  have  been  activated  in  the 
sample  because  the  increase  in  carrier  concentration  at  the 
high  temperatures  gives  a  smooth  curve  rather  than  the 
abrupt  jump  that  is  seen  in  Figure  52.  In  the  one  level 
fit  the  activation  energy  was  0.0537  eV  with  a  donor 
concentration  of  4.76  x  1012  cm--*  and  an  acceptor 
concentration  of  9.13  x  1010  cm'*. 

In  order  to  fit  the  deep  energy  level,  a  two  level  fit 
was  tried.  Although  Chi-squared  was  reduced  to  8.73,  the 
fit  resulted  in  unrealistically  large  values  for  the 
acceptor  and  shallow  donor  energy  level  concentrations. 

The  acceptor  concentration  was  reported  to  be  6.66  x  1013 
cm-3  while  the  shallow  donor  was  7.07  x  1013  cm-3 
with  an  energy  level  of  0.0369  eV.  Because  the  two  level 
fit  was  unsuccessful,  a  three  level  fit  was  tried.  The 
best  fit,  Fig.  54,  gave  a  Chi-squared  of  7.25  with  an 
acceptor  concentration  of  1.32  x  1011  cm-3.  The 
shallowest  energy  level  was  0.0525  eV  with  a  concentration 
4.14  x  1012  cm-3.  The  next  higher  donor  energy  level 
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Figure  54.  Three  Level  Least  Squares  Fit  of  Carrier 
Concentration  of  Si:N  Annealed  at  900  °C 
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was  0.1856  eV  at  1.91  x  10  cm  and  the  deepest  donor 
level  was  0.348  eV  at  6.51  x  1012  cm-3.  Although  the 
0.19  eV  energy  level  was  reported  by  Tokamaru,  this  appears 
to  be  a  coincidence.  An  examination  of  Fig.  54  shows  that 
this  energy  level  just  does  not  fit  the  data.  The  deepest 
level  does  not  fit  the  data  very  well  either. 

The  shallow  level  in  the  three  level  fit;  however,  is 
consistent  with  the  one  level  fit  of  the  FZ0001-A  and  the 
one  level  fit  of  the  unannealed  sample,  FZ0006.  It  appears 
that  this  boule  of  nitrogen  doped  FZ  silicon  has  a  donor 
with  an  energy  level  of  about  0.052  to  0.054  eV  at  a 
concentrat ion  of  4  x  1012  cm-3.  This  energy  level  is 
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not  precise  due  to  the  r-factor  problem;  however,  it  should 
still  be  close  enough  to  identify  the  donor. 


VI .  Analysis 

The  results  of  these  experiments  show  that  nitrogen 
doped  silicon  is  not  as  well  behaved  as  silicon  doped  with 
other  group  V  atoms.  All  of  the  samples  had  some 
temperature  dependence  in  R1/R2  at  low  temperatures. 

Some  samples,  FZ0004-A,  FZ0005,  and  SI-N800C,  were 
temperature  dependent  over  the  entire  temperature  range. 
Because  of  the  requirement  that  R-|  /R2  must  have  a 
constant  value  with  temperature,  it  means  that  these 
samples  are  inhomogeneous.  The  carrier  concentration  data 
from  the  Si:N  samples  that  were  homogenous  enough  for  curve 
fitting  could  only  be  resolved  for  the  shallowest  energy 
level.  The  Si:N  NTD  samples  were  dominated  by  the 
phosphorus  doping  and  this  was  the  only  energy  level  found 
in  the  data. 

The  Hall  experiment  does  not  give  a  precise  activation 
energy  level.  The  r-f actor  problem  distorts  the  energy 
level  even  more.  The  average  donor  energy  level,  0.0483 
eV,  found  in  the  Si:N  NTD  samples  was  close  enough  to 
phosphorus  to  identify  it  as  the  doping  atom  even  though 
the  energy  was  higher  than  the  accepted  value  of  0.045  eV 
(37:21).  The  shallow  energy  level  of  0.052  to  0.054  eV 
found  in  the  Si:N  samples  could  be:  arsenic  which  has  an 
accepted  donor  energy  value  of  0.054  eV,  phosphorus,  or 
nitrogen.  The  donor  atom  is  unlikely  to  be  arsenic  as  no 
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arsenic  was  found  in  the  NTD  samples.  This  leaves 
phosphorus  or  nitrogen  as  the  probable  donor  species. 

Although  Zorin  reported  the  activation  energy  of  ion 
implanted  nitrogen  in  silicon  as  0.045  eV,  no  one  has 
duplicated  this  result.  Zorin's  experiments  were  conducted 
from  100  to  250  K  (44).  Mitchell  repeated  the  Hall 
experiment  on  ion  implanted  nitrogen,  but  extended  the 
temperature  range  from  20  to  250  K  (25).  It  is  probable 
that  Zorin  did  not  find  an  activation  energy  for  nitrogen 
at  all  and  that  his  samples  were  contaminated  with 
phosphorus . 

Mitchell  reported  a  shallow  energy  level  of  0.017  eV 
(25).  Curve  fitting  was  attempted  using  this  energy  level 
on  Si:N  and  Si:N  NTD  samples.  Low  temperature  data  down  to 
20  K  was  included  in  these  attempts  in  spite  of  the 
Rl/R2  problem  without  success.  The  fitting  routine 
kept  pulling  the  energy  up  until  it  closed  in  on  the  0.052 
to  0.054  eV  level. 

A  photoluminescence  (PL)  study  of  unannealed  samples  of 
this  boule  of  Si:N  was  performed  in  1982  by  David  Brown  of 
the  University  of  Dayton  Research  Institute  (5,6).  He 
reported  finding  a  1.1223  eV  line  in  this  material.  This 
line  was  reported  by  Ta j ima  and  is  considered  to  be  an 
exciton  bound  to  substitutional  nitrogen  (38,39).  Figure 
55  is  a  PL  plot  of  the  intensity  verses  photon  energy  of  a 
sample  of  boule  number  0274  (Si:N)  which  was  provided  by 
David  Brown  (6).  The  1.1223  eV  bound  exciton  line  appears 
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Figure  55.  Photoluminescence  Plot  of  Si:N 


between  the  more  prominent  lines  at  1095  and  1135  milli-eV 


which  are  free  excitons  emitted  by  the  silicon.  Brown  also 
found  phosphorus  in  the  samples  at  levels  of  1  to  2  x 
10^2  cm-3.  Other  samples  of  this  material  were  ion 
implanted  with  gallium  and  sent  to  the  MIDAC  Corp.  of  Costa 
Mesa,  California  for  PL  analysis  using  the  processing 
infrared  active  dopant  (PIORAD)  method.  They  also  reported 
finding  phosphorus  and  the  1.1223  eV  nitrogen  line  (6). 

Brown  also  noted  that  although  phosphorus  is  found 
everywhere  in  the  samples,  nitrogen  may  not  be.  When  the 
laser  is  moved  to  another  spot  on  the  same  sample  the 
1.1223  eV  line  may  not  appear  at  all.  This  is  another 
confirmation  that  the  material  is  inhomogeneous. 

Because  the  PL  studies  found  phosphorus  in  the  samples 
in  about  the  same  quantity  as  the  shallow  donor  in  the  Hall 
experiments,  this  positively  identifies  the  shallow  donor 
in  the  Si:N  samples  as  phosphorus.  A  shallow  activation 
energy  level  for  nitrogen  could  not  be  found  in  these  Hall 
experiments.  It  is  possible  that  the  concentration  of 
substitutional  nitrogen  in  bulk  material  is  so  low  that  it 
is  beyond  detection  limits  of  the  system  or  overcompensated 
by  acceptors. 

It  has  been  shown  by  these  experiments,  that  annealing 
activates  deep  levels  in  Si:N.  However,  annealing  may  also 
bring  out  inhomgeneities  that  were  not  present  in  the 
material  before.  The  deep  energy  levels  and  inhomogeneous 
nature  of  the  material  may  be  due  to  silicon-nitrogen 
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complexes  and  precipitates  ( 23 ; 32 : 5020 ) .  Annealing  could 
produce  changes  in  these  complexes.  Ironically,  it  may  be 
the  clustering  of  silicon  nitrides  that  lock  the  lattice 
and  provides  the  strengthening  mechanism. 


VII. 


Conclusions  and  Recommendations 


The  objective  of  this  study  was  to  determine  if  bulk 
grown  nitrogen  doped  silicon  is  a  suitable  VLSI  and  VHSIC 
substrate  material.  These  applications  require  substrates 
with  well  behaved  electrical  parameters.  The  manufacturing 
process  requires  wafers  that  can  stand  up  to  thermal 
processing  without  warping  or  developing  more  defects. 

Small  quantities  of  nitrogen  added  to  the  boule  during 
growth  was  discovered  to  lock  dislocations,  but  the 
electrical  parameters  had  not  been  investigated. 

Findings 

Through  the  series  of  Hall  effect  experiments  performed 
in  this  study  it  was  discovered  that  the  electrical 
parameters  of  this  material  are  not  well  behaved.  The 
experiments  are  listed  in  Table  II  (page  53)  while  Table  IV 
(page  54)  summarizes  the  electrical  parameters  at  300  K. 

The  resistivity  and  the  Hall  mobility  of  the  Si:N  samples 
decreased  upon  annealing.  The  resistivity  in  sample  FZOOOb 
dropped  from  1180  ohm-cm  to  302  ohm-cm  upon  annealing  at 
800  °C  while  the  Hall  mobility  dropped  from  1380 
cm^/v-sec.  Annealing  activated  more  carriers  in  this 
sample.  The  carrier  concentration  went  from  3.85  x  10^ 
em~3  2.97  x  1  0^.  However,  the  decrease  in 
resistivity  and  Hail  mobility  and  increase  in  carrier 
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concentration  is  not  a  function  of  the  annealing 


temperature .  The  sample  annealed  at  900  °C  had  higher 
resistivity  and  Hall  mobility  than  the  samples  annealed  at 
800  °C. 

The  electrical  parameters  in  the  bulk  grown  bi:N  varied 
considerable.  In  all  but  the  unannealed  sample,  the  Hall 
mobility  was  not  as  high  os  the  phosphorus,  doped  CZ 
sample.  The  CZ  sample  had  a  Hall  mobility  et  1300 
cm^/V-sec  while  the  annealed  Si:N  samples  varied  from  881 
crn^/V-sec  to  1  320  cm^/V-sec. 

The  Si:N  NTD  samples  had  similar  resistivity  and 
carrier  concentrations  -if  300  k,  but  not  t  h>  Hill 
coefficient  or  Hall  mobility.  One  Si:N  NTD  sample  had  a 
Hall  mobility  of  1930  cm2/V-sec  and  the  other  sample  was 
1370  cm^/v-sec.  The  high  purity  FZ  NTD  sample  had  a  Hall 
mobility  of  1530  cm^/v-sec  which  is  near  the  mobility 
limit  1500  cm^/v-sec  which  is  due  to  phonon  scattering. 

The  Hall  effect  data  was  taken  over  a  wide  temperature 
range:  20  to  380  K.  Examination  of  the  carrier 

concentration  plots  of  the  phosphorus  doped  CZ  sample,  the 
unannealed  Si:N  sample,  and  the  NTD  samples  shows  that  the 
concentration  decreases  at  high  temperatures.  This 
indicates  that  the  r-factor  or  "r",  is  temperature 
dependent  and  the  assumption  that  it  is  equal  to  one  for 
all  temperatures  is  not  justified.  However,  the  assumption 
that  r  =  1  is  the  best  that  can  be  done  for  now,  because 
the  computation  of  the  temperature  dependent  r-factor  in 
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n-type  material  has  not  been  worked  out  in  detail  (23). 

A  necessary  condition  for  the  accuracy  in  the 
resistivity  and  Hall  coefficient  measurements  is  that  the 
sample  must  be  homogeneous.  The  R1/R2  ratio  was  found 
to  be  temperature  dependent  in  several  samples.  This  means 
that  these  samples  are  inhomogeneous.  Three  out  of  four 
Si:N  samples  were  found  to  be  inhomogeneous  after 
annealing.  The  unannealed  sample,  FZ0006,  became 
inhomogeneous  after  annealing.  The  only  sample  that  was 
homogeneous  after  annealing  was  FZ0001-A  which  was  annealed 
at  the  highest  temperature,  000  °C.  The  Si:N  NTL)  samples 
were  annealed  at  800  °C  but  were  homogeneous.  It  is 
possible  that  these  samples  were  so  dominated  by  the 
phosphorous  doping  that  this  overcame  any  inhomogeneities. 

The  inhomogeneous  nature  of  Si:N  make  it  a  poor  choice 
for  a  VLSI  or  VH5TC  substrate.  The  inhomogeneity  implio:. 
that  the  electrical  parameters  will  not  be  uniform  across 
the  substrate.  This  leads  to  device  problems  such  as 
non-uniform  signal  propagation  time,  different  gains  in 
transistors,  variations  in  resistors,  and  host  of  other 
problems.  The  material  may  be  acceptable  for  discrete  and 
power  devices,  as  these  devices  are  large  enough  to  average 
out  electrical  parameters  across  the  substrate. 

Exactly  what  happens  to  the  nitrogen  in  a  silicon 
lattice  cannot  be  determined  by  the  Hall  effect 
experiments.  The  inhomogeneities  brought  out  by  annealing 
may  be  due  to  silicon-nitride  complexes  and  possible 


precipitates.  It  was  also  concluded  by  examining  the 
concentration  plots  that  annealing  activated  deep  energy 
levels.  This  was  seen  in  all  the  annealed  Si:N  samples 
including;  the  inhomogeneous  samples,  but  not  the 
unannealed  sample,  FZ0006. 

Attempting  to  curve  fit  the  carrier  concentration  data 
to  find  the  donor  activation  energies  and  species 
concentrations  was  clearly  impossible  for  the  inhomogeneous 
samples.  However,  curve  fitting  was  done  on  the  Si:N  NTD 
samples  and  the  two  homogeneous  bulk  Si:N  samples  (FZ0001-A 
and  FZ0006).  The  activation  energy  of  the  deep  level  in 
the  annealed  Si:N  sample,  FZ0001-A,  could  not  be  resolved. 
The  only  shallow  energy  level  that  could  be  resolved  in  any 
of  the  samples  was  phosphorus.  A  shallow  level  that 
unambiguously  identified  nitrogen  was  not  found.  Nitrogen 
was  identified  in  the  Si:N  boule  by  a  photoluminescence 
study  performed  by  David  Brown  (6).  The  conclusion  is 
nitrogen  does  not  go  into  subsitutional  sites  in  silicon  to 
any  appreciable  extent. 

Recommendations 

Based  on  these  experiments  it  appears  that  nitrogen 
doped  silicon  is  not  a  suitable  VLSI  or  VHSIC  substrate 
material  due  to  its  imhomogeneous  nature.  However,  the 
Hall  effect  experiments  on  the  Si:N  NTD  samples  show  that 
it  is  may  be  possible  to  overcome  this  problem. 

It  is  recommended  that  theses  experiments  be  repeated. 


The  Hall  experiment  should  be  performed  on  each  sample 
before  and  after  annealinq  to  qet  a  better  understandinq  of 
the  effects  of  annealinq.  However,  to  determine  exactly 
what  happens  to  nitrogen  in  a  silicon  lattice  a 
Transmission  Electron  Microscope  (TEM)  study  may  be 
undertaken . 

The  experiments  undertaken  in  this  work  gave  data  on 
the  electron  mobility  of  nitrogen  doped  silicon.  The 
mobility  of  holes  in  this  material  can  be  determined  by  ion 
implanting  a  p-type  material  into  Si:N,  the  performing  the 
Hall  effect  experiment.  The  details  of  this  type  of 
experiment  are  given  in  Appendix  A. 

It  is  also  recommended  that  the  temperature  dependent 
r-factor  in  n-type  silicon  material  be  determined  if  n-type 
material  is  to  be  studied  in  the  future.  One  of  the 
results  of  this  would  be  more  accurate  curve  fitting,  thus 
donor  concentrations  and  activation  energies  would  be 
determined  more  accurately. 


Appendix  A 


Hall  Effect  Experiment  with  Bulk  Grown 
Nitrogen  Doped  Silicon  Ion  Implanted  with  Gallium 


In  order  to  study  the  mobility  of  holes  in  n-type 
silicon,  it  is  possible  to  implant  p-type  ions  into  the 
n-type  silicon  lattice.  By  reverse  biasing  the  p-type 
implanted  layer  with  respect  to  the  n-type  substrate,  the 
depletion  region  electrically  isolates  the  ion  implanted 
p-layer  and  the  Hall  effect  experiment  can  be  performed 
(39:336).  As  long  as  the  reverse  bias  leakage  current  is 
less  than  1  or  2%  of  the  source  current,  the  rest  of  the 
current  is  confined  to  the  ion  implanted  p-type  layer  and 
data  can  be  taken  (36:69-72). 

Leakage  and  sample  current  has  to  be  measured  at  each 
temperature.  Leakage  current  increases  with  temperature, 
so  it  is  likely  that  the  temperature  range  of  this  type  of 
experiment  would  be  limited.  Excessive  leakage  current  is 
the  reason  Zorin  had  to  stop  his  Hall  effect  experiments  on 
ion  implanted  nitrogen  in  p-type  silicon  at  250  K  (44:111). 

Wafers  of  phosphorus  doped  CZ  silicon,  Si:N,  and  Si:N 
NTD  were  ion  implanted  with  gallium  at  the  Avionics 
Laboratory  ( AFWAL ) .  Square  samples  were  made  from  the 
wafers.  Four  contacts,  one  in  each  corner,  were  made  on 
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the  ion  implanted  surface.  Contacts  were  formed  by  ion 
implanting  boron  to  degenerately  dope  the  contact  area  with 
a  p-type  material.  Ion  implanting  is  for  making  contacts 
instead  of  using  a  laser  and  spin  on  dopant  because  the 
laser  pulse  penetrates  the  ion  implanted  layer.  In  order 
reverse  bias  the  pn-junction  a  contact  had  to  be  put  on  the 
back  of  the  sample  on  the  n-type  substrate.  Unfortunately, 
this  contact  was  put  on  the  wrong  surface  on  eleven  out  of 
the  twelve  samples. 

Leads  were  soldered  on  the  remaining  sample  which  was 
Si:N  ion  implanted  with  gallium.  The  sample  was  mounted  in 
the  sample  holder  and  tested  at  room  temperature  in  a  Hall 
effect  system  that  had  been  modified.  At  room  temperature 
the  leakage  current  may  be  excessive,  but  the  test  showed 
that  the  leakage  current  was  small  compared  to  the  sample 
current . 

The  system.  Hall  II,  was  modified  to  include  a  power 
supply  to  reverse  bias  the  pn-junction  and  an  extra 
ammeter.  Figure  56  shows  a  schematic  illustration  of  the 
modified  system.  The  high  side  of  the  current  source 
(Keithley  220)  is  connected  to  the  high  side  of  ammeter  A 
(one  channel  of  the  Keithley  619  DVM)  and  the  low  side  of 
the  power  supply.  Ammeter  A  measures  the  current  through 
the  ion  implanted  layer.  Ammeter  B  measures  the  total 
current  through  the  current  source.  The  leakage  current  is 
determined  by  subtracting  current  through  ammeter  A  (sample 
current)  from  the  current  in  ammeter  B  (total  source 
current.  The  ground  on  the  switching  box  is  the  path  for 
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Figure  56.  Schematic  of  the  Hall  Effect  System  Modified 
for  an  Ion  Implanted  Sample  (23) 


the  return  currents.  The  sample  voltage  is  measured  in  the 
same  way  as  described  in  Chapter  IV  for  a  bulk  sample.  The 
voltage  signal  goes  through  high  impedance  unity  gain 
amplifiers.  A,  then  through  the  switching  box  to  the 
digital  voltmeter  (H-P  3455A). 

The  Hall  II  system  is  not  difficult  to  modify  and  the 
room  temperature  test  of  the  Si:N  gallium  implanted  sample 
indicates  that  an  analysis  of  a  p-type  ion  in  Si:N  would  bo 
feasible.  There  are  sufficient  quantities  remaining  of 
gallium  implanted  Si:N  and  Si:N  NTD  to  make  more  samples. 
Because  a  study  of  the  electrical  properties  of  a  p-type 
ion  implanted  layer  in  nitrogen  doped  silicon  has  not  been 
done,  this  could  bo  a  fruitful  thosis  topic. 


Appendix  R 


Graphs  of  Experiment^  FZ0004-A  and  FZ0005 

The  following  figures  are  the  plots  of  the  temperature 
dependent  electrical  parameters  of  FZ0004-A  and  FZ0005.  Tn 
both  cases  these  samples  had  temperature  dependent 
R-)/R2  ratios.  This  means  that  these  samples  are 
inhomogeneous.  The  plots  are  included  to  show  the  results 
of  Hall  effect  experiments  with  inhomogeneous  samples. 
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Figure  57.  Hall  Coefficient  of  FZ0004-A  and  FZ0005 
versus  Inverse  Temperature 


1  9 


2.63  9.40 


6.  6  22- 93  29.69  36-46  43-22 

:00Q/T  T1/K) 


49  -  =>9 


Figure  58.  Resistivity  of  FZ0004-A  and  FZ0005 
versus  inverse  Temperature 
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Carrier  Concentration  of  FZ0004-A  and 
FZ0005  versus  inverse  Temperature 
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Figure  60.  Hall  Mobility  of  FZ0004-A  and  FZ0005 
versus  inverse  Temperature 
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